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<(^ . ABSTRACT 

The ratio of the luminosity of diffuse 511 keV positron annihilation radiation, measured by INTE- 
GRAL in its four years, from a Galactic "positron bulge" (< 1.5 kpc) compared to that of the disk is ^ 
1.4. This ratio is roughly 4 times larger than that expected simply from the stellar bulge to disk ratio 
r^ , of ~ 0.33 of the Galactic supernovae, which are thought to be the principal source of the annihilating 

CIh' positrons through the decay of radionuclei made by explosive nucleosynthesis in the supernovae. This 

JL \ large discrepancy has prompted a search for new sources. 

t*^ , Here, however, we show that the measured 511 keV luminosity ratio can be fully understood in the 

"*T? ' context of a Galactic supernova origin when the differential propagation of these ~ MeV positrons 

j^ I in the various phases of the interstellar medium is taken into consideration, since these relativistic 

positrons must first slow down to energies < 10 eV before they can annihilate. Moreover, without 

propagation, none of the proposed positron sources, new or old, can explain the two basic properties on 

\^ I the Galactic annihilation radiation: the fraction of the annihilation that occurs through positronium 

^T . formation and the ratio of the broad/narrow components of the 511 keV line. 

^i: ' In particular, we show that in the neutral phases of the interstellar medium, which fill most of the 

^^ \ disk (> 3.5 kpc), the cascade of the magnetic turbulence, which scatters the positrons, is damped 

-Y-\ ■ by ion-neutral friction, allowing positrons to stream along magnetic flux tubes. We find that nearly 

1/2 of the positrons produced in the disk escape from it into the halo. On the other hand, we show 

that within the extended, or interstellar, bulge (< 3.5 kpc), essentially all of the positrons are born 

in the hot plasmas which fill that volume. We find that the diffusion mean free path is long enough 

that only a negligible fraction annihilate there and ~ 80% of them escape down into the HII and HI 

envelopes of molecular clouds that lie within 1.5 kpc before they slow down and annihilate, while the 

remaining ~ 20% escape out into the halo and the disk beyond. This propagation accounts for the 

low observed annihilation radiation luminosity of the disk compared to the bulge. 

. , In addition, we show that the primary annihilation sites of the propagating positrons in both the 

S . bulge and the disk are in the warm ionized phases of the interstellar medium. Such annihilation can 

also account for those two basic properties of the emission, the fraction (~ 93 ± 7%) of annihilation 

via positronium and the ratio (~ 0.5) of broad (~ 5.4 keV) to narrow (~ 1.3 keV) components of the 

bulge 511 keV line emission. Moreover, we expect that the bulk of this broad line emission comes 

from the tilted disk region (0.5 < i? < 1.5 kpc) with a very large broad/narrow flux ratio of ^ 6, 

while much of the narrow line emission comes from the inner bulge (i? < 0.5 kpc) with a negligible 

broad/narrow flux ratio. Separate spectral analyses of the 511 keV line emission from these two region 

should be able to test this prediction, and further probe the structure of the interstellar medium. 

Lastly, we show that the asymmetry in the inner disk annihilation line flux, which has been suggested 
as added evidence for new sources, can also be fully understood from positron propagation and the 
asymmetry in the inner spiral arms as viewed from our solar perspective without any additional 
sources. 

Subject headings: (ISM:)cosmic rays — elementary particles — gamma rays: theory — ISM: general 
- nuclear reactions, nucleosynthesis — (stars:) supernovae: general 

1. INTRODUCTION ation at 511 keV from the inner Galaxy has led to ex- 

The discovery (Johnson et al. 1972; Leventhal et al. tensive studies of the possible origin of the annihilating 

1979) of the diffuse Galactic positron annihilation radi- positrons. These studies have shown (eg. Ramaty & Lm- 

' genfelter 1979; Knodlseder et al. 2005) that, of all the 
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source. 

Recent measurements, however, have raised new ques- 
tions about the source. Analyses of the measurements 
by the gamma-ray spectrometer (SPI) on ESA's INTE- 
GRAL spacecraft of the distribution of the diffuse Galac- 
tic positron annihilation radiation have shown that the 
luminosity of a Galactic positron bulge within 1.5 kpc 
of the Galactic Center is ~ 1.4 times that of the Galac- 
tic disk (Knodlseder et al. 2005; Weidenspointner et al. 
2007, 2008a). Their best fit models are either concentric 
shell sources, or Gaussian spheroidal sources, in a bulge 
extending out to ~ 1.5 kpc combined with a thick disk 
with a half thickness of z ~ 0.4 kpc, that extends from 
the Galactic Center to well beyond the Solar distance 
of 8 kpc. They also find a large additional contribution 
from a spherical halo extending beyond the bulge. 

In order to compare with these observations, we use 
these model fitting definitions of the positron bulge, thick 
disk, and halo, following Knodlseder et al. (2005) and 
Weidenspointner et al. (2007, 2008a) in all of the dis- 
cussion that follows, although both have been defined in 
various other ways in the literature. 

1.1. 511 keV Line Bulge/Disk Luminosity Ratio 

The best-fit analysis of the SPI/INTEGRAL first 
year's data (Knodlseder et al. 2005, Shells-|-Dl) gave a 
511 keV line luminosity of Lb - (0.78±0.04) x lO'*^ ^■^q_ 
tons s^^ in a spherical bulge from to 1.5 kpc and an 
old, thick disk luminosity of Ld ~ (0.26 ± 0.07) x lO''^ 
photons s~^, scaled to a Solar distance of i?o ~ 8 kpc. 
This gave a bulge to disk 511 keV luminosity ratio of 
Lb/Ld ~ 3.0 ± 0.8. The subsequent best-fit analyses of 
the first two years' data, which included more coverage of 
the disk (Weidenspointner et al. 2007, for 4 Nested Shells 
and Old Disk), gave a bulge flux of ~ (0.79±0.07) x lO'^ 
photons cm^^ s^^ in spherical shells out to 1.5 kpc, a 
thick disk fiux of - (1.41 ± 0.35) x 10"^ photons cm^^ 
s^^, and a marginally detected spherical halo flux of 
~ (0.86 ± 0.59) X 10^3 photons cm'^ s'^ from 1.5 to 
8 kpc. To infer the corresponding 511 keV line lumi- 
nosities, we use the fiux to luminosity ratios from their 
earlier modeling (Knodlseder ct al. 2005 Table 3), which 
are essentially just 4ttR^ where R ^ Ro for the bulge, 
about 0.57i?o for the thick disk, and about 0.93i?o for 
the halo. From the best-fit fluxes, we infer a 511 keV 
positron bulge luminosity Lb ~ (0.60 ±0.05) x 10^'^ pho- 
tons s^^ in the spherical shells out to 1.5 kpc, and a thick 
disk luminosity of L^, ~ (0.36 ±0.09) x lO'^-'^ photons s^^ 
for & Lb/Ld ^ 1.7 ±0.5. 

Most recently from more than four years of observa- 
tion, Weidenspointner et al. (2008a, model BD) fit the 
511 keV fluxes and luminosities to a two component nar- 
row (FWHM of 3°, or 0.42 kpc) and wide (FWHM of 11°, 
or 1.5 kpc) spheroidal Gaussian distribution and thick 
disk from more extensive observations that increased the 
effective sampling distance to 0.75i?o for the disk. They 
found best-fit fluxes of ~ (0.75 ± 0.09) x 10"^ photons 
cm~^ s^^ from the bulge and ~ (0.94 ± 0.16) x 10^^ 
photons cm^^ s~^ from the disk. Scaled to Ro ^ 
8 kpc, these fluxes give a best-fit 511 keV luminosity 
Lb ~ (0.57 ± 0.08) x lO'*^ photons s~i in the spher- 
ical bulge < 1.5 kpc, and a thick disk luminosity of 
Ld ~ (0.40±0.06) x lO''^ photons s"!. This gives a bulge 
to disk 511 keV luminosity ratio of Lb/Ld ~ 1-4 ± 0.3. 



They also found a best-fit combined halo and bulge lu- 
minosity of ~ (1.55 ± 0.1) X lO**^ photons s~^ for similar 
scaling, and subtracting the bulge luminosity, suggests a 
spherical halo luminosity of Lh ~ 1 x 10**^ photons s~^ 
beyond the bulge, for a total Galactic 511 keV luminosity 
Lg ^ (2.0 ± 0.2) X lO'^^ photons s'^. 

Best-fit SPI/INTEGRAL spectral analyses by Weiden- 
spointner et al. (2008b) also show an asymmetry in the 
positron annihilation line flux from two opposing inner 
disk components at 0° to 50° to either side of the Galac- 
tic bulge. These show an 80 ± 40% excess in the 511 
keV line fiux at the negative longitudes compared to the 
positive. 

Assuming simply that the positron annihilation and 
production rates are in local equilibrium, these obser- 
vations of the luminosities have been taken to imply a 
similar bulge/disk ratio and asymmetric spatial distri- 
bution for Galactic positron production. This suggested 
production ratio and distribution has challenged super- 
nova source models, since it is much larger than that ex- 
pected from the distribution and mean bulge/disk ratio 
of Galactic supernovae (e.g. Knodlseder et al. 2005; Wei- 
denspointner et al. 2008b). As we show, however, these 
spatial properties can be fully explained by positron 
propagation. 

1.2. 511 keV Line Width & Positronium Fraction 

We also consider the other two fundamental spectral 
properties of the Galactic positron annihilation: the ra- 
tio of the broad to narrow 511 keV line emission and the 
fraction (fps) of the annihilation that occurs via positro- 
nium (Ps) formation. As shown by Jean et al. (2006) 
neither of these can be explained by any of the proposed 
positron sources, either old or new, without extensive 
positron propagation, since all of these sources are ex- 
pected to produce most of the positrons in the hot ten- 
uous phases of the interstellar medium, which cannot be 
their primary annihilation site. For this reason alone a 
detailed treatment is required of the production, propa- 
gation and annihilation in each phase and region of the 
Galaxy. 

Studies (e.g. Guessoum, et al. 1991; Guessoum, et al. 
2005) show that positron annihilation in different phases 
of the interstellar medium occurs in differing ratios of di- 
rect annihilation on both free and bound electrons, to in- 
direct annihilation via positronium formation, depending 
on the ionization fraction and temperature. Direct anni- 
hilation produces two photons at 511 keV, while positro- 
nium annihilation produces either two 511 keV photons 
or a three photon continuum depending on the spin state. 
Positronium is formed 25% of the time in the singlet 
state, parapositronium, which annihilates with a mean 
life of ^ 1.25 X 10~^° s into two 511 keV photons, while 
75% of the time it forms in the triplet state, orthopositro- 
nium, which annihilates with a mean life of ~ 1.4 x 10"'' 
s into a three photon continuum between and 511 keV 
(e.g. Guessoum et al. 1991). 

The SPI/INTEGRAL measurements of the ratio of the 
Galactic 27 511 line flux to that of the 37 continuum in 
the bulge and disk shows that the bulk of the positrons 
annihilate via positronium (Ps) with fps ~ 94 ± 6% 
(Churazov et al. 2005), - 92 ± 9% (Weidenspointner 
et al. 2006), and - 95 ± 3% (Jean et al. 2006). The 
weighted mean positronium fraction of ~ 94 ± 4% also 
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allows us to determine the total Galactic positron an- 
nihilation rates, A, in the bulge and disk, since A — 
(e+/75ii)L. From the above we see that {e~^ /j^u) = 
[2(/ps/4) + 2(l-/p^)]^i, which equals 1.69±0.17. Thus 
the inferred best-fit positron bulge and Galactic disk an- 
nihilation rates are Ab ^ (0.96 ± 0.17) x lO''^ e+ s'^ 
and Ad -- (0.68 ± 0.12) x lO'*^ e+ s^^ using the previ- 
ously derived L^ and L/j. In the hot tenuous plasma 
of the halo, however, the positronium fraction depends 
on the refractory grain abundance, ranging from only 
about 18% with narrow (FWHM ~ 2 keV) 511 keV line 
emission, if all the refractories are in grains, to 42% 
with broad (FWHM ~ 11 keV) if the grains all disin- 
tegrated (Jean et al. 2006). Thus we assume a possi- 
ble range of (e''"/75ii) ^ 0.65 ± 0.07 there. From the 
best-fit halo flux, we thus infer a halo annihilation rate 
of Ah ^ (0.65 ± 0.11) x lO^^ e+ s-\ implying a total 
Galactic positron production rate of ^ (2.3 ±0.2) x 10'*'' 
e+ s-i. 

The width of the 511 keV line also depends on the tem- 
perature of the medium and its state of ionization which 
effects the fraction of positronium formed by charge ex- 
change with H, H2, and He by superthermal, ~10 eV, 
positrons as they slow down (Guessoum et al. 1991; 
Guessoum et al. 2005). The prompt annihilation of this 
fast moving positronium produces a characteristic broad 
(FWHM ^ 6 keV) component of the 511 keV line, while 
the subsequent thermal annihilation, either directly on 
free and bound electrons or via positronium formed by 
radiative combination or charge exchange give much nar- 
rower (< 2 keV) lines. Even the annihilation in the hot 
(T ~ lO^K), tenuous ionized medium which would pro- 
duce a broader (^10 keV) line from the hot free elec- 
tions, can instead produce a narrow line if most of the 
positrons annihilate on refractory dust grains (Guessoum 
et al. 1991; Guessoum et al. 2005). The dust, as it 
does with interstellar molecule formation, provides a 2- 
dimensional regime with much higher interaction rates 
for surface chemistry. 

The width of the Galactic bulge 511 keV line mea- 
sured by SPI/INTEGRAL from the first years' data has 
been recently fitted (Churazov et al. 2005; Jean et al. 
2006) by two components with about 67±10% of the 
emission in a narrow line with a width of 1.3±0.4 keV 
and the remaining 33±10% in a broad line with a width 
of 5.4±1.2 keV by positronium formation in flight, giv- 
ing a broad/narrow 511 keV line flux ratio of ~^ 0.5. The 
combined positronium fraction and relative fluxes of the 
broad and narrow 511 keV line emission strongly con- 
strain the bulge annihilation processes, as we discuss in 
detail below. 

1.3. Positron Production & Propagation 

The large difference between the ratio of the 511 keV 
annihilation line bulge/disk luminosity and that of the 
expected positron production by supernovae has led to 
suggestions there is some new, unrecognized source of 
positrons in the Galactic bulge. This seems unlikely, 
however, since various recent reviews (e.g. Dernier & 
Murphy 2001; Knodlseder et al. 2005; Guessoum, Jean & 
Prantzos 2006) of the potential Galactic positron sources 
all conclude that supernovae are still the most plausi- 
ble source and that other suggested sources, including 
cosmic-ray interactions, novae, and various exotic pro- 



cesses, all seem to be weaker and much less certain. 

However, as we have shown (e.g. Guessoum, Ramaty 
& Lingenfelter 1991), the 511 keV annihilation emission 
only illuminates the annihilation sites, not the sources, 
of the positrons, since these positrons must first slow 
down to energies < 10 eV before they can annihilate. 
Thus, we examine in detail the propagation and annihi- 
lation of the /3+-decay positrons in the bulge, disk, and 
halo, and we show here that the expected positron prop- 
agation and annihilation are very different in the bulge 
and disk, and the cloud and inter-cloud environments. 
When these differences are taken into consideration, we 
find that the spatial distribution, the 511 keV line widths 
and the positronium fraction can all in fact be clearly un- 
derstood in the context of a Galactic supernova origin. 

Although positron propagation has previously been 
discussed, the arguments range from the simple sug- 
gestion (Prantzos 2006) that most of the disk positrons 
might diffuse along dipolar field lines and annihilate in 
the bulge, to Monte Carlo transport simulations (Gillard 
et al. 2007) suggesting that they travel less than a few 
hundred parsecs in all but the hot phase. We therefore 
also carefully examine both the physics of MeV electron 
propagation and the observational evidence of such prop- 
agation from extensive measurements of Solar flare and 
Jovian electrons in the heliosphere. 

Here, we examine in detail the mechanisms of the rela- 
tivistic positron propagation in the various phases of the 
interstellar medium. 

Overall, we assume that the Galactic stellar and in- 
terstellar distributions are defined by two superimposed 
systems: a stellar bulge-disk-and-halo system of stars 
and an interstellar bulge-disk-and-halo system of gas and 
plasma, as shown in Figure 1. The stellar bulge and 
disk populations essentially determine the distribution 
of supernovae from the decay of whose radionuclei the 
positrons are born. Following Ferriere et al. (2007), the 
interstellar bulge, disk and halo define the distribution of 
gas and plasma where the positrons die by annihilation. 

As we discuss in more detail later, the stellar bulge 
is confined within the interstellar bulge, which is blown 
out by a bulge wind to about 3.5 kpc in the Galactic 
plane and above the disk it feeds and merges into the 
halo. The stellar disk extends not only throughout the 
interstellar disk out to at least about 15 kpc, but also 
all the way into the interstellar bulge. The interstellar 
disk, however, extends in only to the outer edge of the 
interstellar bulge, and starts in the so-called "molecular 
ring" at around 3.5 kpc, which is defined primarily by 
the two opposing innermost spiral arms. 

Threading the combined stellar and interstellar sys- 
tems, is the Galactic magnetic field, along whose flux 
tubes the positrons travel from their birth to their death. 
We assume that the flux tubes are nearly vertical (e.g. 
Beck 2001), perhaps dipolar, in the inner part of the 
interstellar bulge and are blown out by the bulge wind 
(e.g. Brcgman 1980; Blitz et al. 1993) roughly radially 
in the outer bulge, out to the molecular ring, where they 
begin to be tightly wound into the spiral arms of the 
interstellar disk (e.g. Beck 2001). There, in giant star 
formation regions along these spiral arms, hot, massive 
OB star associations are formed, whose collective super- 
novae generate superbubbles, that blow out plasma as 
well as flux tubes into the overlying halo (e.g. Parker 
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1979). 

Within this framework wc model the various aspects of 
the production, propagation and annihilation of Galactic 
positrons. 

We show that roughly half of the Galactic positrons 
are born in the extended interstellar bulge (e.g. Ferriere, 
Gillard & Jean 2008) in the inner ^ 3.5 kpc of the Galaxy, 
and nearly all result from SNIa occurring in the hot, ten- 
uous high-pressure plasma fills nearly all of that volume. 
We show that ~ MeV positrons created in these SNIa su- 
pernovae are expected to escape from the remnants into 
the surrounding tenuous plasma. There they propagate 
along magnetic flux tubes by 1-dimensional diffusion, res- 
onantly scattered in pitch angle by turbulently generated 
magnetohydrodynamic (MHD) waves that cascade down 
to the thermal electron gyroradius. In these plasmas the 
diffusion mean free path is long enough that only a negli- 
gible fraction annihilate there and ~ 80% of the positrons 
escape down into the dense HII and HI shells of neighbor- 
ing molecular clouds within 1.5 kpc, where they quickly 
slow down and annihilate, while the remaining ^ 20% 
escape out into the halo and disk beyond. 

In the Galactic disk beyond 3.5 kpc, the filling fac- 
tor of hot (~ 10^ K) plasma is more modest (~ 20%). 
Here the plasma is mostly in superbubbles, created by 
the correlated core-collapse supernovae of massive stars. 
Consequently, the bulk of positron production via ^^Al 
decay, synthesized by these stars, occurs preferentially in 
the superbubbles. The larger bubbles blow out into the 
halo, sweeping the magnetic fiux tubes up with them. 
Again the diffusion mean free paths of the positrons are 
such that the positrons escape along the flux tubes, ei- 
ther up into the halo, or down into the warm ionized 
envelopes, surrounding the base of the superbubble gen- 
erated chimneys, where they slow down and annihilate. 

In the warm (^ 10"* K), but essentially neutral HI gas, 
which fills the largest fraction of the disk > 3.5 kpc, we 
show that the cascade of magnetic turbulence, created at 
parsec scales, is damped by friction between the ions and 
neutral gas long before the cascading turbulence reaches 
the positron resonant gyroradii. Thus the positrons suf- 
fer much less scattering, and stream along the flux tubes 
with an isotropic pitch angle distribution. Therefore the 
positrons also mostly stream out into adjoining phases, 
either into the halo or into those warm ionized envelopes, 
separating the ubiquitous warm HI gas from the hot su- 
perbubbles, where they slow down and annihilate. 

Therefore we find that roughly 1/2 of the positrons, 
produced anywhere in the disk >3.5 kpc, slow down and 
annihilate there, primarily in the warm ionized gas, while 
the rest escape from the disk into the halo. This escape 
thus accounts for the low observed annihilation radiation 
luminosity of the disk and explains the observed ratio of 
annihilation rates in the bulge and disk. 

Moreover, the efficient diffusion of positrons out of the 
hot tenuous plasmas in both the bulge and disk, where 
nearly 70% are born, and into the warm ionized envelopes 
and cold neutral interiors of molecular clouds, where they 
annihilate, also explains both the observed positronium 
fraction and broad to narrow 511 keV line flux ratio. 

The positrons escaping into the halo from the bulge 
and disk also lead to the annihilation of roughly 1/3 of 
all the positrons in the halo, even though there in no 
significant production there. 



1.4. Outline 

Although the processes of Galactic positron produc- 
tion, propagation and annihilation are rather straight- 
forward, all the details still make the determination of 
the expected 511 keV line emission both lengthly and 
complicated. 

We consider positron production from the decay of 
three separate long-lived radioisotopes, each synthesized 
by distinct classes of stars: the "''^Ni and ^^Ti, respec- 
tively, by SNIa and SNIp from thermonuclear explosions 
of ^ Gyr-old accreting white dwarves in close binary 
systems, and the ^^Al by young (< 40 Myr), massive 
(> 8 M0) stars, most likely through Wolf Rayet (WR) 
winds and core collapse SNII and SNIbc. The super- 
novae, SNIa and SNIp, from older stellar populations 
are known to occur in galactic bulges and disks, while 
WR stars, SNII and SNIbc, formed in extremely young 
stellar populations, occur in active, or recent, galactic 
star-formation sites, such as spiral arms, molecular-cloud 
complexes, and galactic nuclei. Thus in section 2 we 
briefly re-examine the production of /?^-decay positrons 
by Galactic supernovae and their spatial and temporal 
distribution in the light of recent observations and the- 
ory 

In section 3 we then explore in detail the expected en- 
ergy loss, scattering, and propagation processes of the 
relativistic positrons from these supernovae in the dif- 
ferent phases of the interstellar medium. We consider 
collisionless scattering of these positrons by small-scale 
fluctuations generated by magnetohydrodynamic turbu- 
lence in the ionized interstellar phases, based on exten- 
sive studies and observations of relativistic electrons in 
the analgous turbulence in the interplanetary medium. 
Within this context we adopt a self-consistent transport 
model tied to field-aligned turbulence dependent on the 
properties of the ambient interstellar phases. 

In sections 4, 5, and 6 we determine the expected 
spatial distribution of positron production, propagation, 
slowing down and annihilation in the various phases of 
the positron bulge, disk and halo. In each of these 
we review the fundamental properties of the interstellar 
medium, which constrain positron transport as well as 
leave signatures in the positron annihilation line emis- 
sion features. Here we explore the very large differ- 
ences between the interstellar properties, and the result- 
ing positron propagation and annihilation, in the well- 
studied local interstellar medium compared to those in 
the inner ^1.5 kpc of our Galaxy, where the diffuse 511 
keV emission is most intense. 

Finally, in sections 7 & 8 we give a summary and con- 
clusions, comparing the expected 511 keV line emission 
with the SPI/INTEGRAL observations. We show that 
the observed spatial and spectral features of the Galac- 
tic positron annihilation radiation can all be fully un- 
derstood in the context of a supernova radionuclide ori- 
gin, when the differential propagation of these ^ MeV 
positrons in the various phases of the interstellar medium 
is taken into consideration. 

2. GALACTIC NUCLEOSYNTHETIC POSITRON SOURCES 

Here we briefly re-examine the production of positrons 
by /3^-decay radioactive nuclei produced by different 
types of Galactic supernovae. We also re-examine the 
rates of such supernovae in the Galaxy and the result- 
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ing spatial and temporal distribution of their production 
of positrons in both the Galactic bulge and disk in the 
light of recent observations and theory, in order to deter- 
mine the expected positron bulge/disk ratio of the such 
positron production, Qb/Qd- 

2.1. Positrons from Supernovae 

Positrons resulting from the decay chains of the longer- 
lived radioactive nuclei, '"'^Ni, ^'^Ti, ^^Al, all produced by 
explosive nucleosynthesis in supernovae, have long been 
thought (Colgate 1970; Burger, Stephens & Swanenburg 
1970; Clayton 1973; Ramaty & Lingenfelter 1979; Chan 
& Lingenfelter 1993; Milne, The & Leising 1999) to be 
the major source of Galactic positrons. The total Galac- 
tic positron production rates from these radionuclei de- 
pend on their supernova mass yields M, the supernova 
occurrence rates v and the survival fraction, /, of their 
positrons in the supernova ejecta as it expands into the 
interstellar medium. 



2.1.1. 



'Ni 



The decay chain of •''^Ni^'''®Co^'''®Fe produces a 
positron only 19% of the time, since "''^Ni decays solely by 
electron capture with a mean life of 8.8 days into ^^Co, 
which also decays primarily by electron capture, and just 
19% by P~^ emission into ^^Fe with a mean life of 111.4 
days. Thus the production rate of positrons from "''^Co 

decay is Qge = 0.19/56^56, or 130x 1043/56J^M56 e+ s-\ 
where /56 is the time-integrated survival fraction of ^^Co 
decay positrons in the expanding supernova ejecta, A''56 is 
the production rate of ^^Ni atoms per s, M56 is the ^^Ni 
yield in Mq per SN, and v is the Galactic supernova rate 
in SN per 100 yr. This Galactic rate alone would greatly 
exceed the observed rate of positron annihilation, if all 
of the positrons escaped from the supernova ejecta, but 
most do not. 

The nucleosynthetic yields of ^^Ni, whose decay energy 
powers the visible light from the exploding ejecta, have 
been extensively calculated in all different types of super- 
novae. This is particularly true of the cosmologically im- 
portant Type la thermonuclear supernovae in accreting 
white dwarves, which are best described by the deflagra- 
tion model W7 of Nomoto, Thielemann, & Yoko (1984) 
that gives a ^^Ni yield of 0.58 Mq. The peculiar sub- 
class of Type la(bg), also called Type Ip, supernovae in 
accreting sub-Chandrasakar white dwarves are expected 
to produce 0.44 M0 of ^^Ni from the helium detonation 
model (Woosley, Taam & Weaver 1986). 

Because of the relatively short decay mean life of ^^Co, 
most of these positrons lose their energy and annihilate 
in the ejecta before it becomes either thin enough for 
them to escape or sufficiently rarefied for them to survive. 
Indeed, quantitative studies of the survival of positrons 
in the expanding supernova ejecta (Chan & Lingenfelter 
1993; Milne, The & Leising 1999), for various supernova 
models, have shown that essentially all of the positrons 
produced by its decay in the most frequent, but massive, 
core collapse SNII supernovae slow down and annihilate 
unobserved in their much denser ejecta and do not escape 
into the interstellar medium. 

In particular, Chan & Lingenfelter (1993) have cal- 
culated /56 of positrons from ^^Co decay expected for 
various supernova model distributions of the ejecta den- 
sity, velocity and magnetic fields. They found that in the 



limit where the magnetic field is thoroughly tangled, the 
positron survival fraction for SNIa deflagration model, 
/56, ranged from 0.1% to 2.5% for unmixed versus uni- 
formly mixed ejecta. Alternatively in the limit that the 
magnetic field is fully combed-out (e.g. Colgate, Petschek 
& Kriese 1980) the SNIa /56 was much greater, ranging 
from 5% to 13% for unmixed and uniformly mixed ejecta. 
Following this same procedure Milne, The & Leising 
(1999) have further shown that a mean survival fraction 
/56 ^ 3.5 ± 2% for the deflagration model gave the best 
fits to SNIa light curves at late times (>1 yr), when the 
positrons from ^^Co are the dominant source of energy. 
Using the deflagration model (Nomoto, Thielemann, & 
Yoko 1984) s^Ni yield of 0.58 Mq, we calculate the ex- 
pected positron production rate Q^q = 75 x 10'*'^/56^'/a 



2.1.2. 



44 



Ti 



Chan & Lingenfelter (1993) have shown that the 
positrons from '''^Ti^'^'^Sc^^^Ca decay, which produces 
a positron 95% of the time from /?"'" emission of ^^Sc to 
*'*Ca, is a significant source of the annihilating positrons. 
The Solar system abundance ratio of '*'*Ca to ^^Fe of 
1.23x10^"^ (e.g. Lodders 2003) requires a similar rela- 
tive nucleosynthetic yield ratio of ^^Ti to ^^Ni, since 
these radionuclei are the primary sources of '*'*Ca and 
■"^^Fe (Woosley & Pinto 1988). The much longer ''"^Ti 
decay mean life of 89 yr would allow essentially all (> 
97%) of its positrons to survive in the expanding ejecta 
of even the most massive supernovae (Chan & Lingen- 
felter 1993). Therefore the total Galactic positron pro- 
duction rate from ^^Ti decay can be scaled by the So- 
lar system abundance ratio of '*''Ca to ^^Fe, assuming 
that half of the Galactic ^^Fe is produced by SNIa (e.g. 
Timmcs, Woosley & Weaver 1995) and /44 ~ 1, so that 

O44 ^ 0.95/44^^44, or 1.2 X lO^^i^ia e+ s~'^. 

The most likely sources of these ^^Ti decay positrons 
are the peculiar SNIp supernovae, typified by SN 1991bg 
(Filippenko et al. 1992; Turatto et al. 1996), so that we 
assume the same spatial distribution as SNIa. In particu- 
lar, the primary source of Galactic ^^Ca appears to be He 
detonations in accreting sub-Chandrasakar mass white 
dwarves, whose calculated (Woosley, Taam & Weaver 
1986; Woosley & Weaver 1994) ^''Ti yields of M44 ~ 
O.O2M56 are ~ 16 times the Solar system M44/M56 ratio. 
And such models also give the best fits to the light curves 
of these peculiar supernovae (e.g. Milne, The & Leising 
1999; Blinnikov & Sorokina 2004; The et al. 2006). The 
M44/M56 ratios calculated for other types of supernovae 
are much less than the Solar system ratio. Thus we as- 
sume the spatial distribution of the SNIp is the same as 
that of the SNIa and we also scale the '*'*Ti production 
to that of 5^Ni in SNIa. 



2.1.3. 



26 



Al 



The additional contribution of positrons, produced 
82% of the time from long-lived (r26 ^ 1.04 x 10^ yr) 
decays of ^^ Al to ^^Mg, can be determined much more di- 
rectly from the measured Galactic luminosity of the 1.809 
MeV line, which accompanies that decay. This luminos- 
ity implies a steady-state Galactic mass M26 ~ 2.8±0.8 
Mq of ^^Al (Diehl et al. 2006), or a positron production 
Q26 ~ 0.82M26/26TOpT26 ~ (0.34 ±0.10) x lO''^ e+ s^^. 



6 



Higdon, Lingcnfcltcr and Rothschild 



Thus they can contribute ~ 15±4% of the total Galactic 
positron annihilation rate of ~ (2.3 ± 0.2) x 10^'^ e+ s~^ 
inferred from the best-fit bulge, disk and halo annihila- 
tion luminosity discussed above. 

Knodlscder (1999) has shown that the '^^k\ 1.809 McV 
line luminosity is strongly correlated with the distribu- 
tion of massive stars in the disk, which confirms that it 
is produced in Wolf Rayet winds and/or in core collapse 
supernovae, SNII and SNIb/c. 

With these calculated yields of ^^Ni, '^'^Ti and ^SAl in 
various types of supernovae, we now estimate the total 
Galactic positron production rate from the decay of ra- 
dionuclei synthesized in supernovae in the bulge and disk, 
assuming no significant production in the halo. 



from the decay of supernova produced ^^Co, ^^T^\ and 

26 AL 



Qb+d — Qbe + Q44 

{75f56iyia + l-2l^Ia + 0.3) X 10^3 



Q 



26 



e+s-1 



(1) 



scaled to the mean Galactic occurrence rate, vja of SNIa 
supernovae per 100 yrs. 

2.2. Mean Supernova Rates in Our Galaxy 

Galactic supernova occurrence rates depend on the 
Hubble class of the galaxies and are commonly defined 
(e.g. van den Bergh & Tammann 1991; Cappellaro, Evans 
& Turatto 1999) in units of SNu, equal to 1 SN per 100 yr, 
times a factor of h^L^^i, where ho is the Hubble constant 
in units of 75 km s^^ Mpc^^ and Lq^^ is the blue lumi- 
nosity of the galaxy in units of (10i°L|^)-^ The Hubble 
classification of our Milky Way Galaxy is Sbc (Binncy & 
Merrifield 1998; Kennicutt 2001). From a study of the 
extra-galactic observations Cappellaro, Evans & Turatto 
(1999) find that the SNIa rate is 0.21±0.08 SNu in Sbc 
galaxies. Thus from the Galactic blue luminosity, L^^i 
of (1.9 ± 0.3) X lOi°L0 (Freeman 1985), and this ex- 
tragalactic rate for Sbc galaxies, we estimate the mean 
supernova rate in our Galaxy for SNIa to be 0.40±0.16 
per 100 years. 

However, recent studies of the secular evolution of the 
Galactic bulge suggest that episodic (^Gyr) bursts of 
star formation could cause modest variations in the su- 
pernova bulge/disk ratio. From observations of oxygen- 
rich, cool giants, Sjouwerman et al. (1998) conclude that 
the most recent burst of star formation in the inner bulge 
occurred roughly 1 Gyr ago and produced ~ 10^ SNII 
and SNIb/c supernovae. Such a burst would then be fol- 
lowed by a much more extended period of enhanced SNIa 
and SNIp, occurring with a range of delay times between 
0.1 and 1 Gyr, calculated from evolutionary models by 
Greggio (2005). Thus, ahhough the SNII and SNIbc su- 
pernovae all occurred during the burst of star formation 
a Gyr ago, we would expect the SNIa and SNIp to still be 
occurring at the present time. Using the average ratio of 
SNIa to SNII for our Galaxy of 0.25 (e.g. van den Bergh & 
Tammann 1991; Cappellaro, Evans & Turatto 1999), we 
would expect 2.5x10^ SNIa from this star burst spread 
over ~1 Gyr, or an added contribution of 0.025 SNIa per 
100 years. Thus we estimate the present supernova rate 
in our Galaxy for SNIa to be 0.43±0.16 per 100 years. 

Using this SNIa rate via we then estimate from equa- 
tion (1) the expected mean Galactic positron production 



Qb+d = [(32 ± 13)/56 + (0.5 ± 0.2) 
+ (0.3 ±0.1)] X W^ e+s"^ 



(2) 



Thus, of the Galactic positron production of ^ (2.3 ± 
0.2) X 10^"^ e"*" s~^, inferred from the observed bulge, 
disk and halo annihilation radiation as discussed above, 
26 Al decay accounts for - 0.3 x lO''^ e+ s-\ and '^''Ti 
decay produces ^ 0.5 x 10**^ e+ s^^, nearly all from SNIp, 
leaving ^ (1.5±0.3) x 10**^ e+ s~^, which can come from 
^^Ni decay positrons, if they have a survival fraction in 
SNIa ejecta of /se ~ 5 ± 2%, including the uncertainty in 
the Galactic SNIa rate. A survival fraction of ^ 5%, was 
in fact just what was calculated (Chan & Lingenfelter 
1993 Fig. 3) for the standard deflagration model of SNIa 
(Nomoto et al. 1984 W7) with a combed out magnetic 
field in the ejecta. It is also quite consistent with the 
mean value of ^ 3.5 ± 2% inferred (Milne, The & Leising 
1999) from SNIa light curves at late times, when the 
positrons become a major source of ejecta heating. 

2.3. Spatial Distribution of Galactic Supernovae 

Here we investigate both the stellar bulge/disk ratio 
{B* /D* ^ {vB/vD)ia+ip) of Galactic supernovae and 
their radial and transverse distributions within the stellar 
bulge and disk. Of primary concern are the SNIa and 
SNIp, which appear to be the source of ~ 85% of the 
annihilating positrons. 

We first determine the relative contributions of SNIa 
occurring in the Galactic stellar bulge and disk. It is well 
known that the SNIa occurrence rate is essentially the 
same across the Hubble sequence of galaxy types from 
ellipticals (E) through late-type spirals (Sd), when such 
rates are defined in SNu, the number of supernovae per 
year per unit of blue luminosity of the parent galaxy (e.g. 
Cappellaro, Evans & Turatto 1999). Across this broad 
Hubble sequence, the fraction of galactic blue-band lu- 
minosities contributed by Galactic disks range from zero 
in ellipticals to 95% in late-type Scd galaxies (Simien & 
de Vaucouleurs 1986), yet the SNIa birth rate per unit 
blue luminosity of parent galaxies remains the same. 

Since the SNIa occurrence rates per blue luminosity 
is independent of the relative contributions of galactic 
bulges and disks, we scale the ratio of SNIa birthrates 
in the stellar bulge and disk of our Galaxy by the cor- 
responding ratio of stellar bulge to disk blue luminosi- 
ties, L^ jL^ . In their classic study Simien & dc Vau- 
couleurs (1986) investigated the systematics of bulge- 
to-disk ratios in the blue-band. The mean value of 
their tabulations of L^ jL^ for our Galaxy is 0.25. 
Thus, we expect that the Galactic stellar luminosity 
bulge/disk ratio implies a similar time-averaged, mean 
stellar bulge/disk ratio (j^_b/^'_d)/o ~ 0.25 for Galactic 
SNIa and SNIp. Adding the star burst contribution of 
0.025 SNIa per 100 years in the bulge, which is a 30% 
enhancement over the mean bulge rate of 0.08 SNIa per 
100 years, gives a current stellar Galactic bulge/disk ra- 
tio {vBlvD)ia+ip ^ 0.25 X 1.3 - 0.33 for SNIa and SNIp 
supernovae. 

To model the spatial distribution of SNIa in the 
Galaxy, we use the representations of Dehnen & Binney 



Galactic Positron Annihilation 



7 



(1998) who parameterized the large-scale stellar distri- 
butions within the Galaxy by comparing models to the 
Galactic rotational and other observational constraints. 
Thus the number density of SNIa in the Galactic disk, 
nD*{R, z)ia, and the number density of SNIa in the 
Galactic bulge, nB*{R,z)ia, expressed in cylindrically 
symmetric Galactocentric coordinates, are 



Tin 



_ V 



nB*{R,z)ia^'-^e ^, (3) 



where r] = \J R? + (z/O.G)^, R represents the planar dis- 
tance from the Galactic center in kpc, z represents the 
distance normal to the Galactic plane in kpc, and the 
constant r^ = 1.9 kpc (Dehnen & Binney 1998). Here the 
normalization constants, rioj^ and Uo^ are found from the 
total Galactic SNIa birthrate and the {vBlvD)ia ratio. 
Dehnen and Binney ascertained that the more than 90% 
of the stellar disk population could be represented by Zd 
of 0.180 kpc. They concluded that the most important 
model parameter was the ratio of the disk scale length, 
Rd, to i?0, the distance of the Sun from the Galactic cen- 
ter. They found that the observational constraints were 
best satisfied by any of four models with Rq — 8 kpc 
and i?d/i?o of 0.25, 0.3, 0.35, & 0.4. Here we will use 
their second model with Rd of 2.4 kpc. We also include a 
star burst contribution to the bulge located in the CMZ 
(Sjouwerman et al. 1998). 

2.4. Spatial Distribution of Galactic Positron 
Production 

In Figure 2 we show the relative spatial distributions, 
F{< r), as functions of distance, r, from the Galactic 
center for SNIa residing in the stellar bulge and disk em- 
ploying the relative number densities, equation (3) for 
nB*iR, z)ia and nB*{R, z)ia, expressed a functions of R 
& z in a Galactocentric cylindrical coordinate system. 



FB*i< r) = 271 nB*{R{r), z{r))iar''dr, 
Jo 

Fd*{< r) = 2tt nD*{R{r),z{r))iar^dr, 



FoB*{< r) = 2n noB*{R{r),z{r))ii^ibcr dr. (4) 
Jo 

FoB*{< f) is the relative spatial distribution of mas- 
sive stars, the expected Wolf-Rayet and SNII & SNIbc 
progenitor sources of ^^Al. Although the great majority 
of massive stars, which synthesize ^^Al, are located in 
the outer (> 3.5 kpc) Galactic disk (McKee & Williams 
1997), some massive stars found in the inner (< 0.2 kpc) 
stellar bulge. Since Galactic Lyman continuum radia- 
tion is also produced by a massive star population (e.g. 
McKee & Williams 1997), we use the ratio of Lyman 
continuum radiation emitted by the inner stellar bulge 
to that from the outer disk to estimate the correspond- 
ing ^^Al ratio. Giisten (1989) estimated a luminosity of 
ionizing photons from massive stars in the inner bulge of 
^ 1.5 ± 0.5 X 10^^ photons s~^. Comparing this to the 
estimated (McKee & Williams 1997) total Galactic value 



of ^ 1.9 X 10^'^ photons s""'^, we expect that roughly 10% 
of the massive stars and their generated ^®A1 occur in 
the inner bulge, amounting to ~ 0.03 x 10^'^ e+ s^^. 

For the remaining disk portion we use the spatial distri- 
bution of OB associations, noBiR,z)iiSiibc, of (McKee 
& WiUiams (1997, eq. 35) with Zd = 0.15 kpc to rep- 
resent that of the ^®A1 made by massive stars. These 
spatial distributions have been normalized at large r to 
their expected relative contributions to the total Galac- 
tic positron production rate fractions: Fb* =■ 0.21, Fb* 
= 0.64, and Fqb = 0.15, based on the SNIa and SNIp 
bulge/disk ratio of 0.33 in the 85% of positrons produced 
by them. 

As we see from Figure 2, of the positrons produced 
from SNIa and SNIp in the disk, 13% are born within 
1.5 kpc. Although the best fit disk model used in the 
SPI/INTEGRAL analyses (Weidenspointner et al. 2007, 
2008a) of the annihilation radiation also includes 12% 
occurring within that radius, in order to treat the prop- 
agation of the positrons before the annihilate, we obvi- 
ously need to include them in the total production within 
1.5 kpc. Once we have determined the distribution of 
the subsequent annihilation of all of the positrons, we 
then adjust the estimated bulge (< 1.5 kpc) and disk 
(> 1.5 kpc) annihilation luminosities to compare with 
the SPI/INTEGRAL bulge and disk components, as we 
discuss in section 7. 

Therefore, the positron bulge production (< 1.5 kpc), 
Qb - [Fb* + Fd.{< 1.5) + FoBi< 1.5)]Qb+d - 
[0.21 + 0.64(0.13) + 0.1(0.15)]Qij+D ^ 0.31Qb+d, and 
the disk production, including the OB contribution, 
Qd - [Fd*{^ - FD*i< 1.5)) + 0.9Fob]Qb+d - 
[0.64(0.87) + 0.9(0. 15)]Qb+_d ^ 0.69Qb+d- Thus the 
effective positron production ratio in the bulge and disk 
within and beyond 1.5 kpc is (Qb/Qd) ^ 0.31/0.69 ~ 
0.45. 

From the ratio of positron production rates in the bulge 
Qb and disk Qb wc now determine the differential prop- 
agation of the positrons within and between those regions 
and the halo before they annihilate. In simplest terms 
we define the resulting positron bulge/disk annihilation 
ratio as, Ab/Ad ~ [Pb-.bQb + Pd:bQd)/(Pb:dQb + 
Pd-.dQd), where Pb-.b and Pd-.b, are respectively the 
fractions of the positrons, born in the bulge B and disk 
D beyond 1.5 kpc, which annihilate in the positron bulge 
-B < 1.5 kpc, while Pb-.d, and Pd-.d are the correspond- 
ing fractions of those positrons which annihilate in the 
disk D beyond the positron bulge > 1.5 kpc. As we dis- 
cuss below, however, we need to investigate the positron 
propagation and annihilation in each of the interstellar 
gas phases within each of the regions. 

3. NATURE OF GALACTIC POSITRON PROPAGATION 

We now look in detail at the nature of these positrons, 
and their propagation, slowing down and annihilation 
in the various phases of the interstellar medium in the 
bulge, disk and halo. In particular, we investigate the 
plasma processes that determine the rate of positron 
propagation, drawing upon the extensive studies and ob- 
servations of relativistic electron propagation in the in- 
terplanetary medium. 

3.1. Relativistic (3^ -Decay Positrons from Supernovae 
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Assuming that the dominant source of Galactic 
positrons is the (3~^ decay of ^®Co, ^''Sc and ^^Al from 
supernovae, the detailed studies of their survival in and 
escape from the ejecta show that the great bulk of the 
surviving and escaping positrons are relativistic. 

Positrons emitted in the ^^Co — > ^^Fc decay are dis- 
tributed in kinetic energy with a maximum of 1.459 MeV 
and a mean of 0.630 MeV, and those from "'^Sc -^ ^''Ca 
decay have a very similar spectrum. The expected (Chan 
& Lingenfelter 1993; Milne, The & Leising 1999) mean 
energy of the escaping positrons from the shorter-lived 
(111.4 day mean life) decay of '"'^Co is close to 0.5 MeV, 
reduced from their initial energy by ionization losses in 
the still dense ejecta, while the mean energy of those from 
the longer lived (89 yr ^''Ti dominated meanlife) decay 
of ^''Sc is essentially unchanged at about 0.6 MeV in the 
much less dense ejecta. The mean energy of positrons 
from the decay of the very long-lived (1.04x10^ yr) de- 
cay of ^^Al is likewise about 0.5 MeV. 

In order for these surviving relativistic positrons to an- 
nihilate in the interstellar medium they must first be de- 
celerated to energies, < 10 eV. In order to slow down 
to thermal energies, these ~ 0.5 MeV positrons, must 
run through ^ 0.12 g cm~^ (Berger & Seltzer 1964), or 
ndsd ^ 6 X 10^^ H cm^^, in unionized interstellar gas, 
and this range is reduced by a factor of 1/(1 -f- Xg) if 
the gas is ionized with an ionization fraction, Xg, relying 
solely on coUisional ionization losses. In the process of 
traveling this slowing down distance, dsd ~ 20/ri(l + Xg) 
kpc, the average velocity of these positrons, /3c is ^ 0.7 
of the velocity of light, so their slowing down time, 



tsd ~ dsd/Pc ~ 9 X \Qi'^/n{\ + Xe) yr. 



(5) 



where n is the mean H density in cm^"'. 

Thus we see that since the slowing down distance of the 
radioactive decay positrons in the nominal phases of the 
interstellar medium is > 0.5 kpc in HI clouds with typi- 
cal n ~ 40 cm""^, > 50 kpc in warm media with typical 
n ~ 0.4 cm"'^, and > 3 Mpc in hot phases with typ- 
ical n ~ 3 X 10"'^ cm~^. Since all of these distances 
are far greater than the scales of their corresponding 
interstellar components, these fast positrons, streaming 
along the large-scale interstellar magnetic field, would 
all escape the Galactic disk before they even slow down, 
unless there is sufficient plasma turbulence to efficiently 
isotropize their trajectories. We thus investigate in de- 
tail under what conditions the interstellar turbulence can 
isotropize the trajectories of these relativistic positrons. 

3.2. Nature of the Turbulence which Scatters 
Interstellar Positrons 

In ionized interstellar phases we expect that the in- 
terstellar MeV positrons are scattered by resonant in- 
teractions with ambient turbulent magnetic fluctuations. 
Consequently, the efficiency of these scatterings in turn 
depends on the fundamental properties of the MHD tur- 
bulence at very small spatial scales of the electron gy- 
roradius, r^- However, the basic properties of MHD 
turbulence is poorly understood. In a recent review of 
turbulence and magnetic fields in astrophysical plasmas, 
Schekochihin & Cowley (2007) state "despite over fifty 



years of research and many major advances a satisfac- 
tory theory of MHD turbulence remains elusive. Indeed, 
even the simplest (most idealized) cases are not fully un- 
derstood." 

An important advancement was the realization that in 
many astrophysical sites MHD turbulence is anisotropic. 
Such anisotropy is a well-observed phenomenon in solar- 
wind plasmas (e.g. Belcher & Davis 1971; Matthaeus, 
Golstein, & Roberts 1990; Osman & Horbury 2007) as 
well as in increasingly more accurate, numerical simula- 
tions (Shebalin, Matthaeus, & Montgomery 1983; Maron 
& Goldrcich 2001; Cho, Lazarian, & Vishniac 2002; 
Mason, Cattaneo, & Boldyrev 2006; Perez & Boldyrev 
2008). In a pioneering study Goldreich & Sridhar (1995) 
developed a sophisticated theory of anisotropic incom- 
pressible MHD turbulence. They presented the concept 
of "critical balance" which predicted in inertial regime of 
turbulent MHD flows, first, a filamentary shape for tur- 
bulent fluctuations, via a relationship between the wave 
numbers parallel (||) and transverse (_L) to the direction 

of the local mean magnetic field, kii oc kj^ and, sec- 
ond, the modeled energy spectra of transverse turbulent 

components possesses Kolmogorov (1941) scaling, fc]^ ' , 
observed in low-speed, terrestrial flows. 

At present, a puzzling dichotomy exists between these 
model predictions and a plethora of solar-wind observa- 
tions and numerical simulations. Power spectra of small 
scale electromagnetic fluctuations in solar-wind plasmas 
have been found to have the expected —5/3 spectral in- 
dex (e.g. Horbury et al. 1996; Leamon et al. 1998; Bale et 
al. 2005), but corresponding spectra in many precise nu- 
merical simulations have spectral indices much closer to 
—3/2 (Mason & Goldreich 2001; Beresnyak & Lazarian 
2006; Mininni & Bouquet 2007; Perez & Boldyrev 2008). 
Yet the critical balance prediction about the wavenum- 
ber relationship of turbulent eddies seems to be consis- 
tent with many numerical simulations (Cho & Vishniac 
2000; Maron & Goldreich 2001; Cho, Lazarian, Vish- 
niac 2002; Beresnyak & Lazarian 2006). Unfortunately, 
analyses of single spacecraft measurements are unable to 
determine the full modal three-dimensional wavevector 
spectra of turbulent fluctuations (e.g. Horbury, Forman, 
& Oughton 2005) and thus cannot confirm or refute the 
critical balance prediction about filamentary-shaped fluc- 
tuations (Beresnyak & Lazarian 2008). 

Although much progress has been made in the under- 
standing and modeling of anisotropic MHD turbulence, 
such models, which assume that the frequency of the 
turbulent fluctuations is low compared to the proton gy- 
rofrequency fJp, still have difficulty in representing high- 
frequency coUisionless turbulence evolving on frequen- 
cies, a; « r^e, which seem to be most efficient at res- 
onantly scattering ^ MeV positrons, where Oe repre- 
sents electron gyrofrequency. The application of such 
models of anisotropic turbulence seems to imply that 
coUisionless scattering of MeV electrons would be mod- 
est and MeV electron transport should be dominated 
by streaming at ^ /3c. Yet a great variety of indepen- 
dent investigations of MeV electron transport in the in- 
terplanetary medium irrefutably demonstrate that such 
particles propagate diffusively in the coUisionless solar 
wind plasma, although the nature of the plasma fluctu- 
ations which actually scatter such electrons remains yet 
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unknown. Thus, we employ a simple phenomcnological 
model of collisionless MeV electron scarttering by tur- 
bulent interstellar fluctuations tied heavily to the well 
determined properties of electron scattering in turbulent 
solar- wind plasmas. 

The propagation of MeV electrons has been well stud- 
ied in the interplanetary medium employing measure- 
ments of electrons accelerated by solar particle events as 
well as by the Jovian magnetosphere. The propagation of 
Jovian MeV electrons is, for the most part, dominated by 
diffusion (Ferreira et al. 2001a, b, 2003). Fits to intensity 
and anisotropy time profiles of solar-flare electrons also 
illustrate the dominance of diffusive transport at these 
electron energies, and, consequently, lead to the straight- 
forward determination of scattering mean free paths via 
fits to intensity and anisotropy time profiles (Kallenrode 
1993). 

Bieber, Wanner, & Matthaeus (1996) related such 
phenomenologically-derived scattering mean free paths 
from solar flare particle measurements to the determina- 
tion of simultaneous power spectra of magnetic fluctu- 
ations at MHD spatial scales. Since it is impossible to 
determine the full three-dimensional wavevector spectra 
of turbulent fluctuation from single spacecraft measure- 
ments (e.g. Fredricks & Coroniti 1976), Bieber, Wanner, 
& Matthaeus employed a simple two-component repre- 
sentation of the anisotropic interplanetary MHD turbu- 
lence. They found that their magnetic-field data con- 
sisted of two anisotropic populations, fluctuations with 
large correlation lengths transverse to the direction of 
the mean magnetic field (slab turbulence) and fiuctua- 
tions with large correlation lengths parallel to the direc- 
tion of the mean magnetic field (quasi-two-dimensional 
turbulence). In such slab models the wave vectors of the 
fluctuations are aligned parallel to the direction of the 
mean magnetic field, and a simple interpretation is that 
in the inertial range slab turbulence consists of Alfvcn 
wave propagating along the mean magnetic field. Later 
investigations (Oughton, Dmitruk, & Matthaeus 2004, 
2006; Oughton & Matthaeus 2005) have demonstrated 
that this simple two-component model of anisotropic tur- 
bulence seems to model the basic properties of solar wind 
MHD turbulence. 

Bieber, Wanner, & Matthaeus found that interplane- 
tary turbulence is dominated^ by quasi-two-dimensional 
turbulence, and the mean ratio of the power in slab 
turbulence to the power in two-dimensional turbulence 
is ~0.15. Further, observing that the modeled two- 
dimensional fluctuations do not contribute to particle 
scatterings, they calculated scattering mean free paths 
expected by resonant interactions solely with turbulent 
slab fluctuations via quasi-linear theory (Jokipii 1969). 
They found that such scattering mean free paths agreed 
well with the values determined independently (of the 
analyses of the turbulent magnetic variations) from data 
fits to intensity and anisotropy time profiles. 

^ Recently, Horbury, Forman & Oughton (2005) have investi- 
gated in more detail the properties of anisotropic solar wind tur- 
bulence, implementing a wavelet-based method to track the direc- 
tion of the local mean magnetic field; such an approach reduces 
greatly the noise in the magnetic field measurements. Employ- 
ing the two-component model, they ascertained that their results 
agreed "remarkably well" with the results of Bieber, Wanner, & 
Matthaeus. 



Based on these discussions, we model interstellar MHD 
turbulence as composite slab/two-dimensional turbu- 
lence. We assume here that the ratio of the power of 
interstellar slab turbulence along the mean magnetic field 
to the power of 2-dimensional turbulence normal to the 
direction of the mean field is 0.15 following Bieber et al. 
(1996). Similarly, we assume that the spectrum of in- 
terstellar slab turbulence steepens from a spectral index 
of —5/3 to —3 at a wavenumber kd > ^/rp, where r^ is 
the proton gyroradius, since Leamon et al. (1998) have 
found that at higher wave numbers the spectral index 
of solar-wind variations varied from —2 to —4.4 with an 
average index of —3. 

To quantify such scattering mean free paths for inter- 
stellar MeV positrons we use an analytic approximation 
(Teufel & Schlickciser 2002) calculated for slab MHD tur- 
bulence employing quasi-linear theory. Following Teufel 
& Schlickciser we assume that energy is injected into the 
interstellar medium at an outer spatial scale, lo- For 
klo > 1 and fc < k^, the modeled turbulence possesses a 
spectral index, s. From the discussion above, we assume 
s — 5/3 and at kd = l/fp, we assume the spectral shape 
steepens to s = 3. Teufel & Schlickciser (2002, their eq. 
58), find that for such slab turbulence the electron scat- 
tering mean free path along the magnetic field is 

A|| = I (Bo/SB^^y -r^^(«' ^' *' ^' ^)' (6) 
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where SB^./Sw is the energy density of the turbtflent mag- 
netic fluctuation in the direction of the mean magnetic 
fleld, c is the speed of light, ilg is the gyro-frequency of 
the positron, the mean magnetic field is Bq, Va is the 
Alfven speed, J5o/-\/47rp, p is the ion mass density, re is 
the electron gyroradius; the average positron speed is (3c, 
and K is a dimensionless quantity involving a Gauss hy- 
pergeometric function 2-F1 (Table 3 of Teufel & Schlick- 
ciser 2002). 

Thus, by such one-dimensional (1-D) diffusion the 
positrons in their slowing down time tsd would be 
distributed along a flux tube a mean length Isd ~ 
(2isdA||/3c/3)^/^ in either direction from their point of 
origin. 

3.3. Ion-Neutral Damping 
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The bulk of the Galactic interstellar mass resides in 
primarily neutral phases (HI & H2), concentrated in a 
cloud population (e.g. Tielens 2005). Moreover, a major 
fraction of the disk (> 3.5 kpc) is filled with warm neutral 
HI (Kulkarni & Heiles 1987). Consequently, the nature 
of MHD turbulence in predominantly neutral interstellar 
phases needs to be addressed. 

In partially neutral plasmas, where magnetic forces act 
directly on charges, and indirectly on neutral atoms via 
ion-neutral collisions, the turbulent NHD fluctuations are 
dissipated into heat when ion and neutral velocities differ 
significantly (Higdon 1984; Goldreich & Sridhar 1995; 
Lithwick & Goldreich 2001). Turbulent MHD cascades 
are quenched at scales of roughly the collision mean free 
path of protons with neutral atoms, if the proton collision 
rate with neutral atoms exceeds the eddy cascade rate. 
The mean free path for proton-HI collisions is IpH — 
5 X 10^*^/n cm at 8000 K, or if hydrogen is fully ionized, 
and He is the dominant neutral species, the mean free 
path for proton-He collisions is IpHe ^ 1-5 x 10^'''/n cm 
where n is the H number density, including both ions and 
neutrals (Lithwick & Goldreich 2001). 

Consequently, they showed that MHD turbulence can 
cascade to small spatial scales, much less than these 
damping collision mean free paths, only if the HI frac- 
tion, riHi/n, is less than a critical value. 



riHl/n < fcrit ~ 5{lpH/lo)^^^- 



(7) 



Therefore in predominantly neutral interstellar phases 
turbulent MHD cascades are halted by ion-neutral colli- 
sions at ^ IpH, spatial scales far greater than the gy- 
roradii of MeV positrons. Hence we expect collision- 
less scattering of positrons by such large-scale turbulent 
MHD fluctuations to be very inefficient. Similarly sim- 
ple elastic collisions with ambient electrons, whose mean 
pitch angle scattering is 90° and mean energy loss is 50%, 
have a mean free path much longer than slowing down 
distance due to ionization losses (e.g. Berger & Seltzer 
1964), so they too offer no significant scattering. 

In the absence of such scatterings, ~MeV positrons 
might be expected to generate plasma waves by a reso- 
nant streaming instability, similar to the creation (Kul- 
srud 2005) of ion Alfven waves by streaming relativis- 
tic cosmic ray nuclei. But these do not appear to be 
effective either. These ~ MeV positrons generate res- 
onant whistler waves (Schlickeiser 2002) at wavelengths 
significantly less than the scale of the thermal proton gy- 
roradius, Vp — Up/flp, at frequencies approaching those 
of thermal electrons, fig, where Up is the thermal pro- 
ton speed. The phase speed of whistler waves in the 
direction of the mean field is the electron Alfven speed, 
Vae — {mpl'n^eY Va- If the positron streaming veloc- 
ity, Vs > Vae, and if the whistler waves weren't damped 
so they grew to saturation, then the resulting magnetic 
field fiuctuations in these waves would interact with the 
streaming positrons via quasi-linear wave-particle inter- 
actions that change their pitch angles, reducing Vg to a 
level just less than Vae- However, in predominantly neu- 
tral interstellar phases such whistler waves are also sub- 
ject to severe ion-neutral damping (Kennel 2008 private 
communication) . 



Thus, we expect that positron streaming velocities 
along the magnetic flux tubes are comparable to their 
particle velocities, /3c, similar to the very high drift 
speeds found (Felice & Kulsrud 2001) for cosmic-ray nu- 
clei in warm HI regions. Nonetheless, in turbulent media 
the flux tubes themselves may essentially random walk, 
so that the propagation is described by so-called "com- 
pound diffusion" (Lingcnfcltcr, Ramaty & Fisk 1971). 



4. POSITRON PROPAGATION & ANNIHILATION IN THE 
POSITRON BULGE 

Using these models of relativistic electron transport, 
we now investigate the propagation, slowing down, and 
annihilation of positrons produced in each of the dif- 
ferent phases of the interstellar medium in the positron 
bulge. Galactic disk and halo. The schematic model of 
the bulge, disk and halo is shown in Figure 1. Because 
the positrons, produced by SNIa and SNIp, are widely 
distributed in the bulge and disk (equation 3) , we assume 
that the fraction produced in the various phases of the 
interstellar medium are proportional to their appropri- 
ate filling factors. Although the supernovae themselves 
disturb the local medium, the ~ MeV positrons from the 
decay of radionuclei escape from the remnant into the 
undisturbed surroundings. 

We estimate the relative probabilities, P, of positron 
annihilation within each phase and positron escape into 
neighboring phases by the following procedure. From 
the properties of the medium in each phase, we deter- 
mine the propagation mode, diffusion or streaming, and 
calculate the diffusion mean free path A|| (equation 6) in 
the MHD scattering mode, or the mean velocity, /3c, in 
the unscattered streaming mode. 

Although we have also made Monte Carlo simulations, 
we estimate the propagation of positrons during their 
slowing down in each phase of the interstellar medium, 
using a simple approximation. Since even the average 
properties and structure of the magnetic field especially 
are very poorly known in the various phases, no more so- 
phisticated treatment seems justified. For the 1-D diffu- 
sion in a uniform medium with a mean free path A|| along 
a magnetic flux tube, A^ positrons produced at a point, I 
— and t ~ 0, will be distributed in both directions along 
the flux tube with a density, n(l) = {N /Isdje^'^^^^'"^^^ , 
where Isd = (2A|[/3ctsd/3)"'^/^, by the time tsd that they 
slow down and annihilate. This corresponds to a mean 
positron density n — N/2lsd over the total mean length 
2lsd^ Similarly for positrons streaming at their veloc- 
ity /3c with an isotropic pitch angle distribution they will 
also have a mean density N/2lsd over a total mean length 
2lsd, where Isd ~ Pctsd/2, when they slow down and an- 
nihilate. 

Thus, if the positrons are produced uniformly along 
sonic mean length ^Ib of the magnetic flux tubes thread- 
ing through some phase x, from the above we expect that 
the probability, Px.x, that such positrons will slow down 
and annihilate in that phase before they escape is crudely. 



^x:x ^^ x''B / ^x^sd J^Or x^B ^x ^sd 

and 

Px:x ~ xIb/(xIb +a; hd) for xIb >x hd, (8) 
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where xlsd is the slowing down length in that phase. The 
remaining positrons have a probability, 1 — Px.x, of es- 
caping from phase x into a neighboring phase, and we 
estimate the relative fractions of the escaping positrons 
that go into each of the adjacent phases from simple geo- 
metric arguments, allowing large uncertainties of ± 50%. 

The probability that positrons which diffuse into a 
phase X will slow down and annihilate there is only 
slightly different from those born uniformly within it. 
These positrons are effectively born at their first scatter 
within the phase, roughly within a scattering mean free 
path of the boundary. So if ^Ib <x hd, the probability of 
their slowing down and annihilating in the phase, Px.x is 
the same as for those born uniformly within, ^x ^B/'^xhd- 
But if x^B >x Isd, then Px.x remains at ^ 1/2, since 
after the first scatter near the boundary half would be 
expected to be spread over xhd in both directions along 
the flux tube by the time they slow down and annihilate. 
Thus only about half of them will slow down and anni- 
hilate in phase x, while the other half will escape, or be 
effectively reflected. As we show, however, in the dense 
labyrinth of clouds in the CMZ and tilted disk the half 
of the positrons that diffuse into clouds from the VH and 
HM and arc reflected back out rather than annihilating 
within, quickly diffuse into another cloud where half are 
again reflected, and by the time they have encountered a 
half dozen other clouds only a small fraction, ^ (1/2)^, 
or < 1% remain. As we also show, however, in the diffuse 
disk beyond 1.5 kpc that is not the case. 

From these simple approximations, we estimate the dif- 
ferent propagation fractions P that relate the positron 
production rates Q and annihilation rates A, as defined 
in equation (4) for the expected bulge/disk annihila- 
tion ratio, expanding that equation to explicitly define 
these fractions in the different phases of the interstellar 
medium within the bulge, disk and halo. 

We consider here the Galactic positron bulge as it 
is defined by the spherical component (< 1.5 kpc) of 
the 511 kcV luminosity in the SPI/INTEGRAL analyses 
(Knodlseder et al. 2005; Weidenspointner et al. 2007, 
2008a). Also based on their model fitting we discuss it 
in two parts: the inner bulge region (R < 0.5 kpc) and 
the surrounding outer, or tilted-disk, bulge region (0.5 
kpc < R < 1.5 kpc). 

In order to compare with the SPI/INTEGRAL mea- 
surements, we need to separate the positron produc- 
tion and annihilation in the bulge and disk into five 
regions, shown in Figure 1, defined primarily by their 
galacto-centic radius, although they also have very dif- 
ferent mixes of interstellar phases, as recently reviewed 
by Ferriere, Gillard & Jean (2007). These are 1) the very 
hot, dense inner or nuclear bulge, Bi, (< 0.5 kpc); 2) the 
hot, middle bulge and tilted disk, Bm, (0.5 to 1.5 kpc); 
3) the hot, tenuous outer bulge. Bo between ^ 1.5 kpc 
and the Galactic molecular ring at around ^ 3.5 kpc, 
dominated by the bulge wind (e.g. Bregman 1980; Blitz 
et al. 1993); 4) the predominantly warm neutral outer 
disk beyond 3.5 kpc. Do, and 5) the hot, tenuous halo, 
H, beyond 3.5 kpc and above the disk. The properties of 
each of these regions, which we discuss below, are sum- 
marized in Tables 1 and 2. 

Based on the distributions of the surface density of the 



three positron source components, given in equation (3): 
the bulge and disk distributions of SNIa and SNIp occur- 
rences, and the disk distribution of massive (OB) stars 
(McKee & Williams 1997), we estimated the fraction of 
positron production in each of these regions. The cu- 
mulative production as a function of Galactic radius is 
shown in Figure 2 for each component. We find that 45%, 
38% and 17% of the positrons from the stellar bulge B* 
component of SNIa and SNIp are produced in the regions 
< 0.5 kpc, 0.5 to 1.5 kpc, and 1.5 to 3.5 kpc, respectively, 
plus 10% of the massive star ^^Al component within 0.5 
kpc. Of the positrons from the stellar disk D* compo- 
nent of SNIa and SNIp, 2%, 11% and 30% are produced 
in those same regions, and the remaining 57% are pro- 
duced beyond 3.5 kpc, together with 90% of the massive 
star ^®A1 component. 

Thus the total positron production in the inner bulge 
Bi within 0.5 kpc, the middle bulge Bm between 0.5 and 
1.5 kpc and the outer bulge Bo between 1.5 and 3.5 kpc, 
including both stellar bulge and disk components, and 
the production in the outer disk Do beyond 3.5 kpc, are 

Qs, = [(4.1 ± 1.6)/56 + (0.10 ± 0.05)] X W^ c+s^i 
Qsm = [(5.7 ± 2.3)/56 + (0.09 ± 0.04)] x lO^-'^ c+s-\ 
Qbo ^ [(8.6 ± 3.5)/56 + (0.13 ± 0.05)] x lO'*^ e+s~\ 
Qdo = [(13.7 ± 5.6)/56 + (0.48 ± 0.24)] x lO'^^ e+s-^(9) 



And using the best-fit positron survival fraction, /sg ^ 
5 ± 2% from § 2.2 above, the total positron production 
rates are 



Qs, = (0.31 ± 0.07)] X 10" c ' s 
Qg.^ = (0.37 ± 0.09)] X W^ e+s-i 
Qbo = (0.56 ±0.14)] X lO''^ 
Qdo = (1.16 ±0.20)] X 10"^ t 



e+s-\ 



(10) 



Within each of these regions, y, we determine the ex- 
pected positron production rates in each separate inter- 
stellar gas and plasma phase, x, by the relative filling 
factors of those phases, fyx, times the production rate 
in that region, Qy, such that Qyx ~ fyxQy We then 
estimate the final positron annihilation rates, Ayix', in 
each phase x' within each region y', as the sums over 
x,y,x' and y' of products of the production rates, Qyx, 
times the propagation fractions, Pyx-.y'x', where each is 
the fraction of the positrons born in yx that propagate 
to and annihilate in y'x', as discussed above and given 
for a; : a; in equation (8). 

In the next three sections, we estimate these propaga- 
tion fractions from modelling the propagation, slowing 
down and annihilation of positrons in each of these re- 
gions and phases starting with the Galactic bulge and 
moving outward in the disk and ultimately into the halo, 
which is an important region of propagation and annihi- 
lation, despite its lack of local production. 

4.1. Inner Bulge (R < 0.5 kpc) 

We model the complex interstellar phenomena of the 
inner region of the Galaxy (R < 0.5 kpc), assuming that a 
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very hot (~ 10* K), high-pressure plasma permeates the 
region, pressure equilibrium exists among the interstellar 
phases (Spergel & Blitz 1992; Carral et al. 1994), and the 
properties of the HII, HI, and H2 regions are related to 
each other via the scenario of photodissociation regions 
(PDR) (Tielens & HoUenbach 1985). 

0-star radiation is the primary photo-ionization source 
in the interstellar medium (McKee & Williams 1997). 
Outside of galactic centers, ultraviolet radiation flux at 
an arbitrary location in galactic nuclear regions is gen- 
erated by nearby OB associations, or by single O stars 
distributed randomly throughout the regions (Wolfire, 
Tielens, & HoUenbach 1990). However, in the immediate 
vicinity of a galactic center ionizing radiation is domi- 
nated by the radiation contribution of active galactic nu- 
clei. In the inner Galactic bulge the bulk of the ionizing 
radiation seems to be generated by randomly distributed 
0-stars which photo-ionize the outer layers of nearby in- 
terstellar clouds (Wolfire, Tielens, & HoUenbach 1990; 
Carral et al. 1991). However, in the outer (> 3.5 kpc) 
Galaxy the great majority of O stars are clumped in OB 
associations, which in turn create dense, compact HII re- 
gions, which surround such OB associations (McKee & 
Williams 1997). In the outer galaxy, only, small frac- 
tion (< 15%) of ionizing radiation generated in these OB 
associations escapes to maintain the diffuse HII. 

Thus the photo-ionized, as well as the neutral, phases 
in the centers of galaxies have been modeled as ensembles 
of PDRs (Wolfire, Tielens, & HoUenbach 1990; Carral et 
al. 1991). Each cloud is viewed as having a spherical 
molecular core of radius, r//^ , which contains the bulk 
of the cloud mass at a molecular hydrogen density n/f^ • 
Surrounding the molecular cores are cold atomic HI en- 
velopes, Ar/// thick, which in turn generate infrared con- 
tinua and fine-structure line emissions (Wolfire, Tielens, 
& HoUenbach 1990). Finally, envelopes of photoionizcd 
HII plasma form outer shells, Arnii thick, surrounding 
these HI shells (Carral et al. 1994), ionized by Lyman 
continuum emission from randomly distributed hot, mas- 
sive stars within the central 0.5 kpc. Following Caral et 
al., we assume that these HII envelopes constitute the 
primary photoionized gas component and a schematic 
model of the clouds is shown in Figure 3. Although the 
the hot, massive O stars that generate the HII ionizing 
emission are thought to occur only in the CMZ, we ex- 
pect that clouds in the inner portion of the surrounding 
tilted disk are also likely to be irradiated with sufficient 
flux to maintain HII envelopes. The extent of such irra- 
diation is not known, however, and here we arbitrarily 
assume that it extends only to 0.5 kpc. 

We also assume that the magnetic flux tubes nomi- 
nally pass through the cold clouds in a parallel array, 
but with large turbulent perturbations superimposed, as 
suggested by magnetic flcld in giant molecular clouds. Li 
et al. (2006), using the SPARC 450 fxm polarimetry ob- 
servations, have found that within some such clouds in 
the disk, the direction of the magnetic field is roughly 
correlated with that of the local Galactic field. However, 
extensive Zeeman measurements of molecular clouds also 
indicate strong field perturbations from supersonic mo- 
tions driven by MHD turbulence with energies compara- 
ble to that of the magnetic fields within the clouds (see 
reviews by Crutcher 1999; Falceta-Goncalves, Lazarian 
& Kowal 2008). 



Thus, we assume a nominal field to calculate the base 
mean length of flux tubes through the clouds, and then 
from an estimate of the scale length of turbulent motions, 
we estimate the increase of that mean length resulting 
from the effective random walk or meandering of the flux 
tubes through the turbulent perturbations. The nominal 
mean flux tube length through cold cloud cores is I'g ~ 
4r/3, where r is the core radius, while the nominal mean 
length through an overlying shell of thickness Ar, is I'g ~ 
(4/3)[(r + Ar)^ — r^]/{r + Ar)^. Conservatively taking 
the outer scale of turbulent motion lo as the step size, 
which minimizes the random walk, we expect a mean 
meandering flux tube length, Ib ^ {1'bI^o)1'b- 

In the inner bulge the molecular clouds concentrate 
into two nested disks, the Central Molecular Zone (CMZ) 
and the surrounding Tilted Disk, which extends well be- 
yond 0.5 kpc out to edge of the middle bulge at 1.5 kpc. 
The CMZ is a highly asymmetric region, which is an el- 
lipsoidal disk, ^ 40 pc thick vertically, with lateral axes 
of ^ 500 pc by ^ 250 pc and an interstellar H2 mass 
of 1.9x10^ M0(Ferriere, Gillard, & Jean 2007). In this 
region we estimate a mean molecular cloud core mass, 
Mn^core ^ 5 X 10^ M0 and a mean core radius, rH2 ^ 5 
pc, from integrations of the mass spectrum (2x 10^ M© < 
M//2core < 2x10^ M0) and the size spectrum (3.3 pc 
(resolution limit) < r^j^ < 10^ pc) of the CMZ molecu- 
lar cloud population determined by Miyazaki & Tsuboi 
(2000). This implies a total of ^ 400 such cloud cores 
with a H2 density n/^j ~ 1900 cm^"^. It also gives a total 
volume filling factor of these H2 clouds in the CMZ is 
.fH2 ^ 0.08, but in the fifll inner bulge < 0.5 kpc volume, 
the filling factor is only //fj ^ 4 x 10^''. 

The much larger Tilted Disk, which is described in 
more detail in the next subsection (4.2), surrounds the 
CMZ, extending out to ^ 1.5 kpc with an estimated in- 
tersteUar H2 mass of 3.4x10^ Mq (Ferriere, Gillard, & 
Jean 2007). About 20% of the Tilted Disk lies within 
the bulge source region < 0.5 kpc. This inner portion of 
the Titled Disk contains ~ 7 x 10^ M© of interstellar H2. 
We assume that this H2 gas is also concentrated in sim- 
ilar molecular cloud cores, which would number ~ 140, 
if they have the same mean density and mass as those 
in the CMZ, and they would have an even smaller filling 
factor in the bulge source volume < 0.5 kpc. 

Within and surrounding the CMZ, containing the 
molecular clouds, is a very hot tenuous plasma (VH). 
Observations of 6.7 keV line emission from the Ka He- 
like transition of ionized Fe indicate that this VH plasma 
has a temperature of ~ 10* K, and fills most of this 
central region (Koyama 1989; Yamauchi et al. 1990; 
Spergel & Blitz 1992; Koyama et al. 1996; Muno et 
al. 2004). Thus we assume that a hot (~ 9 x 10^ K), 
plasma (n ^ 0.04 cm""^) fills the inner ~ 500 pc of the 
bulge (e.g. Yamauchi et al. 1990). The thermal pressure, 
Pth/k — nvHTvH/lJ' ^ 10^ K cm^'^, where /i represents 
the mean particle mass for fully ionized H, and k is Boltz- 
mann's constant. This value agrees well with previous 
estimates of the CMZ thermal pressure of ~^ 5 x 10^ K 
cm^3 by Spergel & Blitz (1992). 

Analyzing the diffuse nonthermal radio emission of the 
Galactic center region, Spergel & Blitz (1992) found that 
the magnetic field pressure is in approximate equilibrium 
with the thermal gas pressure. Consequently, we assume 
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that both the intercloud and cloud field strengths are 
^ 100/iG, the value expected from such pressure argu- 
ments, 5x10^^" erg cm^*^. 

Between the VH plasma and the H2 core of each cloud 
are the two surrounding envelopes, the upper warm, ion- 
ized HII gas and the lower, cold neutral HI gas. The 
HII density is estimated from pressure balance between 
HII gas and the VH (Carral et al. 1994). Although the 
thermal pressure of the overlying VH is approximately 
matched by the total pressure of the cold medium (H2), 
the latter pressure is dominated by turbulence (Tielens & 
HoUcnbach 1985; Spergel & Bhtz 1992; Oka et al. 1998). 
Similarly, the large line widths observed in the extended 
HII gases in the Galactic center region (e.g. Rodriguez- 
Fernandez & Martin-Pintado 2005) imply that turbu- 
lence also dominates the HII pressure. We assume that 
the magnetic field strengths are the same in the VH 
plasma and the HII gas. Thus we find that a turbulent 

HII pressure, pturb = \phii{u^hii) (Tielens & Hollen- 

bach 1985), constituting ~ 85% of the total HII pressure 

equal to the total VH pressure, possesses a root mean 

1/2 

square turbulent velocity, (u'jjjj) , of > 25 km s^^, for 

a mass density, pnii, corresponding to uhii ~ 100 cm^"^. 
Note that the HII thermal pressure, 2nHiiTHii, ^ 10® 
K cm"'^, is small compared to that of the VH, since Thii 
is - 5000 K (Mezger & Pauls 1979). 

Based on Wolfirc et al. (2003), we expect that at the 
high thermal pressures found in the Galactic nuclear re- 
gion, warm HI is unstable and only cold HI is thermally 
stable. From analyses of fine structure emissions from 
PDRs of CMZ molecular clouds located far from thermal 
radio sources Rodriguez-Fernandez et al. (2004) found 
typical HI densities, uhi, of ~ 1000 cm""^ with Thi ^ 
150 K. For R < 0.5 kpc, Ferriere, Gillard, & Jean es- 
timated an interstellar HI mass of 1.7x10® M©, with a 
space averaged density, < uhi >^ 25 cm^'^ in the CMZ 
and a filling factor there of fni ^ 0.025. At such a mean 
density, the typical thickness of cold atomic shells is ex- 
pected to be Arni ^ 0.5 pc. The typical thickness of 
HII shells, surrounding such HI shells with a filling factor 
Jhii ~ 0.045, is expected to be Avhjj ~ 0.7 pc. Thus, 
the outer radius of a typical, double-shelled CMZ cloud, 
Tc ~ 6.2 pc. The total cloud filling factor in the CMZ, 
fc = .fH2 + Ihii + Ihi is ^ 0.15 and it further implies a 

— 1/3 

mean distance between cloud centers dec ~ 2rc/c of 
only ~ 23 pc and that between cloud surfaces, dec — 2rc 
of just ~ 11 pc. 

As we will show, the annihilation of the positrons takes 
place almost entirely in the cloud shells. However, only 
a negligible fraction (< 1%) of them are produced there, 
since the filling factors of the cloud phases are all neg- 
ligible compared to the inner bulge volume. Thus we 
assume that all of the inner bulge positrons are born in 
the pervasive, very hot medium, Qbivh ~ Qsi ~ (0.31± 

0.07) X 10^3 e+ S-l, QB^HU ~ QB^HI - Q B^H, ~ 0. 

The subsequent positron annihilation in these phases 
of the bulge is then defined by the propagation fractions, 
which we estimate below using this general model of the 
CMZ. Because, as we show, the positron slowing down 
and annihilation in the HII phase of CMZ is so efficient, 
the probability of positrons, formed in the VH, pene- 
trating through the HII into the underlying HI and H2 is 



negligible, so we only show here the terms for diffusion 
between the other adjacent phases with significant filling 
factors. We also find that the probability of positrons 
successfully escaping out of the inner disk is quite small. 
Thus the significant terms in each phase are, 

AsiVH ~ PvH-.VHiQsiVH + Q BmPBm-.BiV H 

+ QBoPBo:BiVH), 

ABiHII ^ PHII:HIl{QBiHII + QBiVHPvH-.HII 

+ QBmPBm:BiHII + Q BoPBo-.BiHIl) , 

A-BiHI ^ ABiH2 ^ 0, (11) 

where, like those discussed above, the propagation frac- 
tions P, subscripted VH:VH is the fraction of the 
positrons produced in the VH phase that slow down and 
annihilate in that phase, VH:HII is the fraction that es- 
caped from it into the adjacent HII, and other fractions 
are similar. The propagation fractions and the resultant 
annihilation rates are listed in Table 2. 

4.1.1. Very Hot Medium (VH) 

The bulk of the positrons from the inner bulge SNIa 
and SNIp are expected to be born in the VH, which in- 
cludes both the high density and temperature plasma in 
the CMZ and all of the lower density and temperature 
hot plasma that fills nearly ^ 100% of the inner bulge < 
0.5 kpc. To determine the propagation fractions within 
the VH, and between it and the HII regions, we estimate 
the median distance from the positron production site to 
the HII envelopes around clouds in the CMZ and inner 
Tilted Disk. We consider separately the two distributed 
source components: the compact 57% of the inner bulge 
SNIa and SNIp, resulting from the star burst and in- 
ner disk population, and the diffuse 43%, resulting from 
the supernovae in the general bulge population that are 
expected to occur throughout the inner bulge. 

For the star burst and disk component with positron 
production uniformly distributed sources in the CMZ 
with a mean distance between HII shells would he d ^ 
10 pc, suggests that in a fairly regular magnetic field the 
mean flux tube length Ib would also be ~ 10 pc in the 
CMZ. On the other hand, for the bulge component with 
the production roughly uniform throughout the spheri- 
cal volume, the positrons would be distributed along a 
mean distance d ~ (4/3)7rr3/27rr2 ~ (2/3)r ~ 300 pc 
both above and below the CMZ with similar mean flux 
tube lengths, Ib, along the roughly vertical magnetic flcld 
in that region (e.g. Beck 2001). 

In order to determine what fraction of the positrons, 
born in the VH, can slow down and annihilate there and 
what fraction escape into the adjacent HII envelopes, we 
also need to estimate the diffusion mean free path of these 
positrons in the VH plasma. 

It is difficult to quantify the properties of hypotheti- 
cal turbulent flows in the VH, since the VH origin and 
age are unknown. However, a plausible constraint on 
such turbulence is that the dissipation rate of the tur- 
bulent energy must be less than the cooling rate of the 
VH plasma. The rate, at which turbulent energy is dis- 
sipated is ^ pvHiuyj^Y^'^/lo, in units of erg cm""^ s~^, 
where pvH is the mass density of the VH, (uyj^Y^^ is 
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the root mean square turbulent velocity, and lo is the 
outer scale of turbulence (Townsend 1976). 

From the X-ray observations Muno et al. (2004) esti- 
mated that the VH cooling time, tcooi ~ 10^ yr. Since the 
VH thermal energy density, eth ~ fpt/u the cooling con- 
straint becomes (uyfj)^^'^ ^ {loe-th/ {tcooipVH)Y^^ ■ Tak- 
ing the maximum estimate of lo ~ 10 pc, the mean 
separation between cloud surfaces, this relation limits 
(uy^)^/^ to less than 50 km s"^, which is much less than 
the VH Alfven speed, Va ~ 920 km s~^. This cooling re- 
lation places a severe constraint on VH turbulence, if the 
VH is long lived on the time scale of tcooi- However, if 
the VH is a transient phenomenon, this constraint could 
be weakened and the turbulence would be stronger. 

Thus, if the mechanism generating the VH turbulence 
operated over a duration, t„i, less than tcooi, the dis- 
sipation rate of the turbulent energy can be greater. 
Thus (^)V2 



{lo€thl{t^PVH)f'^- If, 



25 Myr, then (uy^)^/^ 



IS 



~ Eth/U. 

very conservatively, i, 

less than 90 km s~^. Since in strong turbulent MHD 

flows, where energy is equipartitioncd between magnetic 

and velocity fields < SBy^^ >^ AirpvH < Uy^j > s 

(e.g. Iroshnikov 1964; Kraichnan 1965; Biskamp 2003), 

constraining < Uy^ > limits < SBy^ >, and then 

< SBf >i/2, since < 5BJ >^'^^ 0.15 < SB^jj >^'\ 

Thus {Bo/5B\\)^ ^ 800, and, consequently a rather large 



A 



25 pc in the CMZ from equation (6). 



In the VH the positron slowing down time is tgd ^ 

9 X 10'*/n(l + Xe) ~ 1.1 X 10^ yr from equation (5) at n ^ 
0.04 cm^^. In the slowing down time, tgd, the positrons 
diffusing along a flux tube in the VH will be distributed 
over a mean length Isd ~ (2A||/3cisd/3)^/^ ~ 2000 pc, 
since the average positron velocity during deceleration is 
~ 0.7c. 

Roughly 61% of the inner bulge positrons are expected 
to occur in the CMZ, predominantly from the SNIa and 
SNIp in the star burst and the inner disk, along with 
10% of the massive OB stars producing ^®A1. The re- 
maining 39%, all from SNIa and SNIp, are expected to 
occur throughout the nuclear bulge < 0.5 kpc. This mean 
diffusion length is much larger than the mean flux tube 
length, Ib, for positrons escaping into the HII. For those 
produced by the star burst and disk component of SNIa 
and SNIp in the CMZ and inner disk, where the mean dis- 
tance between the HII envelopes of neighboring clouds, 
Ib is only ^ 10 pc, we would expect from equation (8) 
that PvH-.VH ~ Is/'^lsd < 0.01, so only a negligible frac- 
tion would slow down and annihilate in the VH, while 
nearly all would diffuse into the HII. 

The much more extended positron production by the 
bulge supernovac fills the inner bulge more uniformly 
within 0.5 kpc. Within the full volume the mean density 
of the VH is closer to ^ 0.007 cm"'^ and the temperature 
is much lower, ~ 6 x 10^ K (Almy et al. 2000). Assum- 
ing a mean Bo ~ 17/xG in equilibrium with the thermal 
pressure, and taking the same tm of 25 Myr, and a longer 

10 ~ 50 pc, we expect a mean free path A|| ^ 30 pc from 
equation (6). Thus, with a slowing down time from equa- 
tion (5) oftsd ~ 6x 10^ yr, the slowing down length Igd ^ 
5000 pc, while the median flux tube length of Ib ^ 300 
pc between the CMZ and the middle bulge, so Pvh-.vh ^ 



0.03. 

Thus, weighting these two probabilities with their rel- 
ative (0.61:0.39 including ^^Al) production fractions, we 
would expect the overall Pvh-.vh ~ 0.02. Monte Carlo 
simulations give the same result. They also suggest that 
less than 2% of the positrons escape beyond the bulge 
(> 1.5 kpc), since they tend to be reflected, because 
the scattering mean free path there is much shorter and 
the slowing down time is much longer, so there is little 
likelihood of their slowing down and annihilating there 
before they diffuse back out and eventually diffuse into 
the HII below. Although, as we discussed above, about 
half of the positrons entering the turbulent HII envelopes 
are also reflected, because of the closeness of clouds, 
positrons escaping from one cloud envelope quickly en- 
ter another, so after just half a dozen cloud encounters 
nearly all have slowed down and annihilated in the HII. 
Therefore, wc conclude that the effective propagation 
factors are Pvhvh ^ 0.02 and Pvhhii ^ 0.96, while 

PvH-.Be - 0.02. 

Thus, only 2% of the positrons born in the VH and 
diffusing along meandering magnetic flux tubes, would 
be expected to slow down and annihilate there before 
they moved into the HII shells surrounding the molecular 
clouds. 

4.1.2. Photoionized Medium (HII) 

From the source- weighted filling factors we expect that 
only ~ 0.03 of the SNIa and SNIp positrons are born in 
the warm ionized HII phase, but over an order of magni- 
tude more positrons are expected to diffuse into it from 
the surrounding VH. 

The highly turbulent nature of the HII medium follows 
from observations of the CMZ that reveal extreme radio- 
wave scattering by small-scale, turbulent density fluctua- 
tions (Lazzio & Cordes 1998; Goldreich & Sridhar 2006). 
Although the origin of such intense small-scale turbu- 
lence is still uncertain, such small-scale fluctuations are 
expected to occur in the ionized interfaces between the 
dense, molecular cloud cores and the VH (Lazio & Cordes 
1998; Goldreich & Sridhar 2006). 

As discussed previously, the HII shells around molec- 
ular cloud cores have a density nnii ~ 100 cm""^ and a 
thickness, Arnii ~ 0.7 pc. The electron fraction high, 
Xe ~ 1, in view of the expected presence of the small- 
scale, ionized density fluctuations. Thus, a positron 
slowing down time, tsd, is ~ 450 yr from equation 
(5). We model the spatial scale of the HII turbulence, 
lo ^ 0.25ArHii ^ 0.2 pc, assuming that the properties 
of such turbulence can be approximated by those gener- 
ated by shear flows in channels Arnii wide (Pope 2000). 
Further, wc find that (u^jj)^^^ ^ 25 km s^^ from our as- 
sumption that turbulent pressure is about 6 times that of 
the thermal pressure. From equation (6) we expect that 
A|| is very small, ~ 1.7 x 10~^ pc. In the slowing down 
time, tsd, these positrons diffusing along a flux tube in 
the HII plasma will be distributed over a mean length 
;,d - (2A|t/3ct,d/3)i/2 - 0.3 pc. 

This positron slowing down distance is less than the 
expected nominal mean length of a flux tube through 
the shells, I'b ~ (4/3) [(r + Arf - r^]/ir + Arf - 1.2 
pc, and much less than that expected for the meandering 
flux tubes in the turbulent ionized gas, which effectively 
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random walk on the scale of the turbulence, Ig- For a tur- 
bulent scale equal to a quarter of the HII shell thickness 
or ^ 0.2 pc, as discussed above, the mean meandering 
flux tube length would be Ib ~ Q-'b/^o)1'b) ~ '''•2 pc, or 
more than 20 times the slowing down length. 

From the filling factor we expect only a negligible frac- 
tion of the positrons to be born in the HII, while, as we 
saw, essentially all (^ 96%) of the positrons born in the 
inner bulge, are expected to diffuse into the HII envelopes 
from the VH. There in each encounter, as we discussed 
above, roughly 50% are expected to quickly slow down 
and annihilate within the outer Isdi which is much less 
than the mean flux tube length Ib through the HII. Since 
these positrons are effectively born close to the boundary, 
the other half are expected to escape back into the VH. 
But there, because of the close proximity of the clouds, 
the escaping positrons quickly diffuse into another HII 
envelope, where the process is repeated, and soon virtu- 
ally all are slowed down and annihilated in the HII. Thus 
the net propagation fractions are effectively Prii-.hii ^ 
1 and Phii-.vh ^ Phii-.hi ^ 0. 

Thus we expect essentially all of the positrons that are 
either horn in the inner bulge or diffuse into it from be- 
yond, slow down and annihilate in the HII envelopes, 
making them the annihilation trap of the region. 

4.2. Middle Bulge (0.5 < R < 1.5 kpc) 

As Ferriere, Gillard, & Jean (2007) have pointed out, 
"reliable observational information" in this region of the 
Galaxy is very sketchy and the gas distributions are very 
uncertain. The middle bulge from 0.5 to 1.5 kpc is dom- 
inated by a tilted elliptical interstellar gas disk, observed 
primarily in CO and HI lines, with a hole in the middle 
of sufficient size to envelop the CMZ (Ferriere, Gillard, 
& Jean 2007) . They suggest that the most plausible disk 
model is that of Liszt & Burton (1980) with a semi-major 
axis of 1600 pc and a semi-minor axis of 600 pc, and a 
central hole are 800 pc and 260 pc, all tilted by 29° rel- 
ative to the Galactic plane. Thus, as can be seen in 
our Figure 1, this disk appears in a more face-on pro- 
jection on the plane of the sky and extends ^ 800 pc 
above and below that plane. This is twice the nomi- 
nal 400 pc scale height of the overall Galactic thick disk 
modeled in the SPI/INTEGRAL analyses of the 511 keV 
line emission (Weidenpointner et al. 2007, 2008). Thus, 
we might expect that much of emission from positron 
annihilation in the Tilted Disk would not have been 
counted as part of the modeled thick planar disk fitted 
from the SPI/INTEGRAL data and would instead have 
been counted as part of their modeled extended spherical 
bulge from which it is essentially indistinguishable. 

Here we model the interstellar phenomena of this outer 
portion of the tilted interstellar gas disk, again assuming 
that a hot tenuous plasma permeates the region, pressure 
equilibrium exists among all the interstellar phases, mag- 
netic pressure is in equilibrium with gas thermal pressure, 
and the properties of the HI, and H2 regions are related 
to each other via the PDR scenario. 

This disk contains an estimated ~ 3.4 x 10^ M0 of 
H2 and the corresponding disk-averaged gas density is 
< nH2 >~ 4.8 cm"3 (Ferriere, Gillard, & Jean 2007). 
In the 80% of the disk beyond 0.5 kpc, we assume that 
the H2 also resides in compact molecular clouds which 
are roughly similar to those in the CMZ with a nominal 



molecular core radius of 5 pc, but because of the lower 
external pressure, these clouds have a lower mean H2 
density n//^ ~ 1000 cni^^. Thus, the typical core mass, 
MH2core IS ^ 2.6 X 10^ M0, and the number of such 
cloud cores, NH2core, is ^ 1300. For this H2 component 
we expect a cloud H2 filling factor in the outer part of 
the tilted disk, fn^ ~ 4.8/1000 ~ 0.005. 

PDR shells surrounding molecular cloud cores must be 
sufficiently thick (~1.5 A„) to shield the underlying H2 
gas from dissociation via FUV absorption (Wolfirc, Hol- 
lenbach & Tielens 1993). A visual extinction, A^, ~ 1.5 
corresponds to an HI column density N{HI) ^ 3 x 10^^ 
cm~2 (Tielens 2005). Thus to shield the tilted disk 
molecular cloud cores from dissociation, uhi^hi of the 
overlying HI shells, ~ 3 x 10^^ cm^^. If we assume that 
the nominal value of the HI density uhi ^ 1000 cm^^, 
half the density of the H2 cores, then Arni ^ 1 pc. Thus, 
for the HI interstellar component we expect a filling fac- 
tor, fni ~ [(1 + ArHi/rn^f - 1] x (4.8/1000) - 0.004 
in the disk. 

These clouds, like those in the CMZ, are likely to be 
highly turbulent with the turbulent pressure balanced by 
the thermal pressure of the external hot medium (HM). 
Also as with the CMZ, we rely on analyses of X-ray ob- 
servations to supply us with the expected HM pressure. 
In their analysis of the ROSAT 0.75 kcV all-sky survey, 
Almy et al. (2000) found that this region of the Galactic 
bulge is filled with a hot (< 10^ K), tenuous (< lO'^ 
cm~^) plasma. However, this thermal HM pressure is 
significantly greater than the thermal pressure at which 
both cold HI and warm HI can coexist, based on Wolfire 
et al. (2003). Thus, we expect that the cold HI is the 
only stable HI phase in the tilted disk region 

Unlike the CMZ and inner tilted disk region, the outer 
tilted disk beyond around 0.5 kpc does not appear to 
have any significant HII component. Although a popu- 
lation of molecular clouds exists in the Tilted Disk, no 
giant HII regions, standard tracers of massive star forma- 
tion, have been resolved in this region (Smith, Biermann, 
& Mezger 1978). Moreover, no stars younger than 200 
Myr have been found here by van Loon et al. (2003) in 
the mid-IR survey data from DENIS and ISOGAL. So 
we conclude that the stellar Lyman continuum flux and 
its resulting HII component are negligible in this region. 

Therefore the filling factors in the tilted disk are ///^ ~ 
0.5%, fni ~ 0.4%, and fnM ~ 99%. The net cloud filhng 
factor (/c = fni + ///2 ) implies a mean distance between 

neutral cloud centers dec ^ '2rc/fc of ^ 60 pc, and a 
mean distance from a random point in the disk intercloud 
medium to a cloud surface of {dcc/2) — {re + ArHi) ^ 25 
pc. 

We do not expect the SNIa and SNIp to be concen- 
trated in this highly tilted disk which extends up to ^ 
800 pc above the Galactic plane. For as we saw from 
equation (3) and Figure 2, 55% of these supernovae are 
in the stellar bulge population, distributed throughout 
the spherical shell from 0.5 to 1.5 kpc, and the other 45% 
are in the disk population, concentrated in the Galactic 
plane with a scale height of ~ 0.18 kpc. Thus, the median 
distance from positron sources in the HM to the labyrinth 
of HI cloud envelopes in the disk is ~ 0.5 kpc for both 
populations. With the roughly vertical magnetic field in 
that region (e.g. Beck 2001), the positrons are born along 
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flux tubes above and below the disk with a mean length 
/s ~ 1 kpc, twice the median value. We also expected 
that essentially all of the positrons are produced in the 
HM with the filling factor in the entire middle bulge, 
Ihm ^ 1, and the other filling factors in the full middle 
bulge volume are fH2 ^ 5 x 10"^ and fni ^ 4 x 10^'"'. 
Thus Qhm -- Qsm -^ (0.37 ± 0.09) x iO'*^ ^+ g-i^ ^nd 
Qh2 ^ Qhi ^ 0. 

Since the HM filling factor is essentially unity, the sub- 
sequent positron propagation and annihilation in each of 
the phases of the outer titled disk is determined predom- 
inantly by the propagation fractions, which we estimate 
below, for the diffusion of positrons from the HM into 
the HI shells and on into H2 cloud cores, where they 
slow down and annihilate. 

Ahm ~ QhmPhm-.hmi 
Ah I ^ QhmPhm-.hiPhi-.hi, 

Ah2 ^ QhmPhM:HiPhI:H2Ph2:H2i (12) 

where, like those discussed above, the propagation frac- 
tions P, subscripted HM:HM, HLHI and H2:H2, are the 
fractions of the positron that slow down and annihilate 
in that phase, independent of whether they are produced 
in that phase or escaped into it from an adjacent phase, 
as given by the other P combinations of escape from one 
phase to another adjacent phase. These fractions and 
annihilation rates are also listed in Table 2. 

4.2.1. Hot Medium (HM) 

Since the hot medium fills effectively all of the middle 
bulge, we expect that all of the positrons produced be- 
tween 0.5 and 1.5 kpc are born in the HM. Based on 
the analysis of Almy et al. (2000) at a Galactic pla- 
nar radius of ^ 1 kpc, which is the midpoint of this 
region, hot gas temperature Thm ~ 5 x 10® K, the den- 
sity uhm ~ 5 X 10~^ cm~^. Thus the positron slowing 
down time tsd ~ 9 x 10® yr from equation (5). 

As with the VH, we estimate the diffusion mean free 
path of the positrons in the HM as follows. We assume, 
as did Spergel & Blitz (1992) for the CMZ, that the mag- 
netic field pressure is in equilibrium with the gas pres- 
sure. Thus we expect Bo ^ y/Sirpth ^ 12 fiG. As with 
the CMZ a plausible constraint on the nature of such tur- 
bulent flows is that the rate at which turbulent energy is 
dissipated is less than the rate at which the hot plasma 
cools. Almy et al. (2000) estimate the cooling time of 
this hot plasma, ic ^ 3 x 10* yr, but we again consider 
a conservative value of ^ 5 x 10^ yr. Approximating an 
outer scale, lo ^ 50 pc for a hypothetical turbulent HM 
flow, this cooling time limits (u^)^/^ to be less than 45 km 
s""'^ for a thermal energy density, eth ~ l-5pt/i ~ 10~^^ 
erg cm"'^. Since the estimated Alfven speed, Va ^ 310 
km s~^, this limit on (u^)^/^ corresponds to A|| ^ 45 
pc by equation (6). In the slowing down time, tsd of 
9 X 10® yr, HM positrons can diffuse a mean distance 
hd ~ (2A||Md/3)i/2 ~ 7 kpc. 

This distance is much larger than the mean flux tube 
length in the HM, /s ~ 1 kpc, along which the positrons 
are born above the labyrinth of HI covered clouds in 
the tilted disk. Thus as we saw with the VH, we ex- 
pect from equation (8) that the fraction of positrons 



born in the HM and slowed down and annihilated there, 
Phm-.hm ~ Is/^lsd ^ 0.07. As Monte Carlo simula- 
tions also suggested with the VH, only a few percent of 
the positrons are likely to escape into the tenuous region 
above 1.5 kpc, while the bulk (~ 91%) diffuse into the 
cloudy labyrinths of both the tilted disk and the CMZ be- 
low. There due to the closeness of nearby clouds within 
these disks, essentially all of the positrons from the HM 
that enter these two labyrinths rapidly diffuse on into 
the turbulent outer envelopes of clouds where they slow 
down and annihilate. 

We further estimate that the relative fractions of 
positrons born in the middle bulge that enter into the 
labyrinth of clouds in the outer tilted disk and annihi- 
late in their HI envelopes, compared to those that escape 
down into the inner bulge and enter into the labyrinth of 
clouds in the CMZ and inner tilted disk and annihilate 
in their HII envelopes, are proportional to the relative 
surface areas of the two disks, or about 75% to 25%, re- 
spectively. Thus, of the 91% of the positrons from the 
HM that enter these regions, we expect Phm-.hii ~ 0.22, 
diffusing into the inner bulge, and the ~ 0.69% remaining 
in the middle bulge are all expected to diffuse into the 
HI shells of clouds in the titled disk, and as we discuss 
below, essentially all slow down and annihilate in those 
shells, giving an effective Prmhi ~ 0.69 and Phmh2 ^ 
0. 

4.2.2. Cold Medium, (HI) 

As with the molecular clouds in the CMZ, we assume 
that only a fraction of the pressure is thermal in the 
tilted disk molecular cloud cores, so the mean temper- 
ature Th2 ~ 20 K and density nj/j ^ 1000 H2 cm^^, 
although we assume they also have nominal radii of ~ 5 
pc. But because of the negligible filling factors of these 
cloud cores, no significant positron production should oc- 
cur within them, and we show that essentially all of the 
positrons born elsewhere are expected to slow down and 
annihilate in the overlying HI envelopes before they can 
get into the H2 cores. 

In the overlying HI envelopes between the cold H2 
cloud cores and the hot HM plasma, we assume a mean 
density of nni ^ 1000 H cm^^ at a temperature Thi ~ 
100 K, giving it thickness Ar/// ~ 1 pc, as discussed 
above. Therefore the positron slowing down time in 
the HI, tad is ~ 90 yr from equation (5), and since the 
MHD cascade is also damped by ion-neutral friction, the 
isotropic streaming positrons are distributed along a flux 
tube over a mean length I ad ~ Pet ad /'^ ^ 9 pc in both 
directions in their slowing down time. 

This is longer than the assumed nominal mean length 
(2/3) [(r + /\rf - r^]/{r + Ar)^ - 1.7 pc of un- 
perturbed parallel flux tubes passing through a HI shell. 
But it is less than the mean meandering flux tube length 
^B ^ {1'b/^o)1'b) "^ 12 pc, expected for an turbulent scale 
lo ^ 0.25ArHi ^ 0.25 pc, if we assume the same shear 
flow scaling (Pope 2000), as discussed above. 

Only a negligible fraction (^ 0.4%) of middle bulge 
positrons are expected to be born in the HI envelopes, 
but as we saw, ~ 69% of all these positrons are expected 
to diffuse into the HI from the HM. As with the positrons 
that diffuse into the HII envelopes of CMZ clouds, these 
positrons are all effectively born within a scattering mean 
free path of the surface. Thus, since lad < ^s, we expect 
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also that ^ 0.5 slow down an annihilate in the HI, be- 
fore they can penetrate into the H2 cores below, and the 
remainder escape back into the HM above. There, again 
because of the close proximity of the clouds, the escaping 
positrons quickly diffuse into another HI envelope, where 
the process is repeated, and soon virtually all are slowed 
down and annihilated in the HI. Thus the net propaga- 
tion fractions are effectively Phihi ^ 1 and Privh ^ 
0. 

Although the density of the HI envelopes are quite un- 
certain, the expectation that they are thick enough to 
slow down and annihilate the positrons is rather robust. 
That results from the fact that the thickness of the HI 
and the slowing down length are both inversely propor- 
tional to the density, because the column depth of the 
HI, which shields the H2 from dissociation, is fixed and 
so is the required column depth for the slowing down. 
Therefore, Z^, lad, lo, Ib and Is/hd all nominally scale 
proportionally. 

Thus, we expect most (^ 70%) of the positrons that are 
either born in the moddle bulge or diffuse into it from the 
outer bulge, slow down and annihilate in the HI envelopes 
of the clouds in the tilted disk, making them the annihi- 
lation trap of the region, as the HII envelopes are in the 
inner bulge. 

As we discuss in the summary (§7), nearly all of the 
positrons that slow down in the HI envelopes in the mid- 
dle bulge annihilate via fast positronium formation in 
charge exchange, producing broad 511 keV line emission, 
while those that slow down in the HII envelopes in the 
inner bulge annihilate predominantly via thermal radia- 
tive positronium formation, producing narrow 511 keV 
line emission. 

4.3. Outer Bulge (1.5 < R < 3.5 kpc) 

Based on the stellar disk distribution of SNIa and 
SNIp, as discussed above, we see from Figure 2 that ^ 
30% of all the stellar disk positrons and ~ 22% of all 
the stellar bulge positrons produced by them should be 
made in the outer bulge between 1.5 kpc and 3.5 kpc, 
and, as we now show, they become a major contributor 
to the modeled positron bulge component of the observed 
annihilation radiation within 1.5 kpc. 

This region has been known for some time (Faber & 
Gallagher 1976; Giisten 1989; Dehnen & Binney 1998) to 
contain little HI, as traced by 21 cm radio line emission 
(Lockman 1984), and little H2, as traced by 2.6 mm line 
of CO (Sanders, Solomon & Scoville 1984; Robinson et 
al. 1988). The lack of neutral gas results from a "high 
pressure galactic wind" (e.g. Blitz et al. 2007, Breg- 
man 1980), driven primarily by stellar bulge supernovae 
within ~ 0.5 kpc, that blows out to at least 3 kpc, before 
it can recombine to form new neutral gas. 

In this region, therefore, we consider only the hot 
plasma, with a mean density n ^ 2 x 10^'^ cm^"^ av- 
eraged over the scale height of ^ 2 kpc, a conservative 
mean magnetic field strength, _Bo ^ 4 fiG, an outer scale 
for MHD turbulence, Iq ^ 75 pc, based on simulations of 
supernova-driven MHD turbulence (de Avillcz and Bre- 
itschwerdt 2004), and, finally, we assume that the root 
mean square strength of the turbulent magnetic field 
SB /Bo ^ 0.5. From equation (6), we then find that the 
scattering mean free path parallel to the mean magnetic 
field. All ^ 14 pc. Assuming that the bulge wind blows 



the magnetic field out radially, the mean length of a fiux 
tube through the region is Ib ^ 2 kpc. 

In this ionized plasma the slowing down time, tsd ^ 
23 Myr from equation (5). In that time the positrons 
diffuse both directions along a flux tube a mean length 
Isd ~ (2A||/3ctsd/3)^/^ ^ 6 kpc. Comparing this length 
with that of the fiux tube, we expect from equation (8) 
that the probability that positrons, born uniformly in 
the Bo, slow down and annihilate there is ^ Is/'^lsd ^ 
0.17. But since the contribution from this region was 
also included in the model halo of Weidenspointner et al. 
(2008a), we include Bo there instead and take 1.5 kpc as 
the inner boundary of the halo. 

The remainder escape both into the bulge below and 
the disk and halo above. As we see from the distribu- 
tion of the supernovae in Figure 2, the 70% of the Bo 
supernovae in the disk are rather uniformly distributed 
between the two boundaries and we assume that equal 
fractions escape up and down. However, since the re- 
maining 30% from the bulge component are concentrated 
along an Zb < 1 kpc close the inner boundary at 1.5 kpc, 
we expect that they escape primarily into the bulge and 
have a smaller chance of annihilating in the Bo. Thus, 
we estimate that the net fraction of Bo positrons anni- 
hilating there is Pbo-.Bo ~ 0.14 and that escaping down 
into the middle bulge PBo-.Bm ~ 0.57. 

As with the positrons produced in the middle bulge, 
we further expect that the ^ 57% of the outer bulge 
positrons going down into the Bm from the Bo will also 
meet the same fates in the same proportions, mostly dif- 
fusing into the labyrinths of the tilted disk and the CMZ, 
annihilating in the HI and HII envelopes of the molecular 
clouds. Thus, of all the positrons produced in the Be we 
expect Pbo-.hi ^ 0.40 and Pbo-.hm ~ 0.04 in the middle 
bulge and Pbo-.hii ~ 0.13 in the inner bulge. 

Those positrons going beyond 3.5 kpc are also expected 
to annihilate and not return, and within large uncertain- 
ties we assume that they have roughly equal likelihood 
of going into the halo or the disk. Thus Pbo-.Do ~ 0.14, 
while since the Be anniilation was also included in the 
halo, the net Pbo-.h ^ 0.29. 

4.4. Positron Bulge (< 1.5 kpc) Annihilation 

We now estimate the total positron annihilation rate 
within the 1.5 kpc bulge, using the contributions of the 
various phases and the propagation fractions within and 
between them, summarized in Table 2. 

Using these fractions and combining equations (9) 
through (13), the total annihilation rate in the positron 
bulge < 1.5 kpc and its components are thus expected to 
be, 



Ahii ^ (0.46 ± 0.08) X 10'" e ' s 
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Ahi -- (0.48 ±0.17) X 10'" e's 

AvH+HM -^ 0.06 X 10^3 e+s-i 

AB^ -- (0.47 ± 0.08) X lO''^ e+s"i 

ABm -- (0.53 ±0.13) X 10''^ e+s"i 

Ab ~ (1.00 ±0.15) X 10 



« e+s-i. 



(13) 



As we discuss in detail in the summary (§ 7), this is 
quite consistent with the best-fit bulge component anni- 
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hilation rate from the SPI/INTEGRAL analyses (Wei- 
denspointner et al. 2008a). Although we assumed large 
uncertainties of ±50% in the estimated propagation frac- 
tions between interstellar gas phases, which produce 
comparably large uncertainties in the individual anni- 
hilation rates in each, their sums are more tightly con- 
strained. 

We also see (Table 2) that although the hot phase of 
the interstellar medium fills most of the bulge, because 
of propagation it is not the site of most of the annihi- 
lation. We find (Table 2 and equation 13) from the ex- 
pected propagation and slowing down, that ^ 46±8% of 
the bulge positrons annihilate in the HII shells surround- 
ing the molecular clouds in the highly ionizing radiation 
environment of the inner bulge < 0.5 kpc, via thermal 
positronium formation resulting in the narrow (^ 1.2 
kcV), 511 keV line emission (Guessoum, Jean & Gillard 

2005) from the bulge. We also see that ~ 48±12% of the 
annihilation is expected to occur in the middle bulge, 
0.5 to 1.5 kpc, in the neutral HI shells surrounding the 
molecular clouds which lack significant HII shells, via in- 
flight positronium formation resulting in the broad (~ 
5.8 keV), 511 keV line emission from the bulge. As we 
discuss in detail below (§7), these two annihilation sites 
can account for the observed broad/narrow 511 keV line 
ratio of ~ 0.5 (Churazov et al. 2005; Jean et al. 2006) in 
the positron bulge. We further point out that since these 
two components are expected to be essentially formed in 
separate, resolvable regions, spectral analyses of the two 
regions can test this prediction and further probe the 
structure of the bulge interstellar medium. 

Only a fraction Pbo-.Bo ^ 14% of the positrons arc 
expected to slow down in the outer bulge, producing 
an annihilation rate Abo ~ 0.09 x 10**^ e+ s""'^ in an 
essentially spherical shell contribution between 1.5 and 
3.5 kpc. This is consistent with the effective limit of 
< (0.2±0.2) X 10**^ e"*" s~^ from that particular region, in- 
ferred from the best-fit modehng of the SPI/INTEGRAL 
data (Knodlseder et al. 2005). But since the contribution 
from that region has subsequently been included as part 
of the model halo in the analyses of Weidenspointncr et 
al. (2008a), as discussed above, we include Abo m the 
halo annihilation instead. Even though the annihilation 
fraction in the outer bulge is small, its contribution to 
the observed 511 keV line flux is still significant, because 
of the high 75ii/e+ yield of ~ 1.7 expected (Jean et al. 

2006) from positron annihilation in dusty hot plasma, as 
discussed above (§ 1.2). 

In the next section we follow the remaining ~ (0.16 ± 
0.10) X 10^'^ e+ s"^ of the outer bulge positrons that 
escape beyond 3.5 kpc into both the halo and the inner 
spiral arms that deflnc the Galactic molecular ring of the 
interstellar disk. As we also show below (§ 5.4), the inner 
disk annihilation dominated by the these arms, which 
are viewed asymmetrically from our solar perspective, 
can fully account for the disk flux asymmetries (Weiden- 
spointncr et al. 2008b). 

5. POSITRON PROPAGATION & ANNIHILATION IN THE 
GALACTIC DISK (> 3.5 KPC) 

Having estimated the expected production, propaga- 
tion and annihilation of positrons in the positron bulge 
(< 1.5 kpc) and in the outer bulge from 1.5 to 3.5 kpc, wc 
now turn to the various components of the positron pro- 



duction, propagation and annihilation in the outer disk 
{Do) beyond 3.5 kpc. 

As we have seen, the total expected positron produc- 
tion rate in this region consists of those born there, QDo^ 
and those that diffuse into it from the outer bulge be- 
tween 1.5 and 3.5 kpc, QboPbo-.Do, we estimate 



Qdo-- (1.16 ±0.20) X lO'^^ e+s-\ 
QboPbo-.Do - (0.08 ± 0.04) X 10^3 ^+g-i_ 



(14) 



We again focus on the major phases of the interstellar 
medium, which fill essentially all of the outer disk vol- 
ume in which essentially all of the supernovac in this re- 
gion occur, namely: the warm predominantly neutral HI 
medium, (HI), the hot tenuous medium (HT), and the 
photoionized medium produced by O stars: the dense 
dense HII regions (HII) and the diffuse warm ionized 
medium (WI). The dense ionized phases together with 
the molecular clouds, OB star associations and super- 
bubbles are all concentrated along the spiral arms (e.g 
Drimmel & Spergel 2001; Stark & Lee 2006), which dom- 
inate the structure of the outer disk. Their properties, 
discussed below, are summarized in Table 3. 

We assume nominal average values of filling factors 
of these phases beyond 3.5 kpc of fni ~ 0.50, /ht ~ 
0.20, Jhii ~ 0.15 and fwi ^ 0.15, based roughly on 
recent models (Ferriere 1998, 2001) of a thick disk of ^ 
400 pc half thickness. There are large uncertainties in 
these numbers, but since they provide only an estimate 
of where the positrons are born, the uncertainties do not 
have a significant impact on where they eventually slow 
down and annihilate, for as we have seen elsewhere, that 
is dominated by their propagation. Although we do not 
consider independently the other interstellar phases of 
the disk, namely the molecular H2 and HI clouds, be- 
cause of their very small filling factors (< 1% Ferriere 
1998), we do consider their effects in the diffuse neutral 
medium (HI). 

As we saw from the integration of equation (3) and 
Figure 2, 57% of the SNIa and SNIp in the disk occur > 
3.5 kpc, producing ~ (0.86 ± 0.30) x 10^^ g+ g-i^ ^g 
assume that they occur randomly in the various phases 
of the disk proportional to their filling fractions and that, 
as discussed above, - 90% of the - (0.3 ± 0.1) x 10^^ e+ 
s^^, produced by ^^Al decay from massive stars in the 
disk, are born in the HT phase of hot superbubbles. Thus 
we expect that the positron production rates in each of 
the major phases of the disk beyond 3.5 kpc are 



Qht-- (0.45 ±0.12)] X 10 
(0.45 ±0.15)] X 10'*-'' 



4-^ c+s-i, 



iHI 



c ' s 



Q 



HII 



}wi ^ (0.13 ±0.04)] X 10*^ e 



(15) 



From the positron production in these phases, we es- 
timate the propagation fractions, listed in Table 4, and 
the resulting positron annihilation in the disk beyond 3.5 
kpc as follows. 

Art ^ Pht:Ht{Qht + QhiiPhii-.ht + QwiPwi-.ht), 
Ami ~ Phi:Hi{Qhi + QhiiPhii-.hi + QwiPwi-.hi), 
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^Hii ^ Phii:Hii{Qhii + QhiPhi-.hii + QwiPwi-.Hii 

+QhtPht:Hii + QboPbo-.hii, ) 

Awi ~ Pwi:Wi{Qwi + QhiPhi-.wi + QhiiPhii-.wi 

+ QHTPHT:Wl)il&) 

where, like those discussed above, the propagation frac- 
tions P, subscripted HT:HT, HLHI, HILHII, and WI:WI, 
are the fractions of the positrons that are either produced 
in that phase or escape into it from an adjacent phase, 
that then slow down in and annihilate in that phase, 
while the other Ps are all combinations of escape from 
one phase to another adjacent phase. That excludes only 
direct escape from the HT into the HI, or vise versa, be- 
cause they are separated from one another by the HII 
and WI, and any return from the halo, which we assume 
is negligible. 

Because of the large vertical scale height of the hot 
tenuous plasma in the large blown out superbubbles, an- 
nihilation there is also effectively part of the halo, not 
a resolvable component of the ±400 pc disk, so we add 
Aht to Ah- 

5.1. Hot Tenuous Medium (HT) 

Ferriere (2001) found that the mean HT filling fac- 
tor in the outer disk was a modest, Jht ^ 0.2. She 
also found that the HT was dominated by the contri- 
bution of superbubbles (e.g. Tomisaka & Ikeuchi 1986; 
Heiles 1987, 1990; Mac Low & McCray 1988) created by 
core-collapse supernovae, which in turn are correlated in 
space and time. We further found (Higdon & Lingen- 
felter 2005) that the majority (^ 85%) of core-collapse 
supernovae (SNII & SNIb/c) occur in superbubbles and 
an additional fraction (^ 20%) of the SNIa and SNIp 
occur there randomly with a probability equal to the su- 
perbubblc filling factor. 

Positron production in the superbubbles therefore 
comes both from the decay of ^® Al from the massive star 
Wolf Rayet winds and SNII, which create them, and from 
the decay of ^^Co and ^^Ti from the random SNIa and 
SNIp. Thus, as discussed above, the production in the 
HT is expected to be Qht ~ (0.45 ±0.12) x lO^^ g+ g-i^ 
with Q26 ^ (0.27 ± 0.09) x lO^^ 3+ g-i f^om clustered 
massive stars and Qia+ip ~ (0.18 ± 0.04) x 10'*'^ c+ s^^ 
from the random SNIa and SNIp. 

The typical properties of the hot medium (e.g. Ferriere 
1998; Yan & Lazarian 2004) are: a nominal tempera- 
ture Thai ^^ 1 x 10^ K, a scale height averaged density, 
nHM ^ 2 X 10""^ cm"'^, and a mean magnetic field, Bq ^ 
4 /xG. As discussed in §4.1.1, the rate at which turbu- 
lence is dissipated must be less than the radiative rate 
at which the plasma cools. However, the cooling time 
of HM in outer Galaxy is very uncertain. Therefore we 
used tcooi ^ 3 X 10^ yr, the log average coohng time for 
a hot medium, whose evolution is dominated by time- 
dependent cloud evaporation following the prescription 
of McKee & Ostriker (1977), and a cooling time from 
Tielens (2005) for the above values of Thm and uhm- 
Using an outer scale, lo ~ 75 pc (de Avillez and Bre- 
itschwerdt 2004), the above cooling time limits (u2)i/2 



to be less than 40 km s'^ With V„ ^ 80 km 



,-1 



and 



the root-mean-square SB / Bq < 0.5, equation (6) implies 
A|| > 10 pc. 



In this very tenuous plasma the positron slowing down 
time is tsd ^ 23 Myr from equation (5) and the slow- 
ing down length Igd ^ {2X\\f3ctsd/3Y^^ > 6 kpc. As- 
suming a typical flux tube length of /^ < 1 kpc in 
the superbubbles blowing out into the halo, the effec- 
tive propagation fraction, Pht-.ht ^ Is/^lsd < 0.08. 
Since, as noted above, these scales are much larger than 
the disk supernova distributions, there is no resolvable 
disk component of annihilation in the HT. Thus the 
Aht ^ Pht-.htQht < 0.04 x 10*^ e+ s"^ is included 
in the halo Ah- 

As with the bulge wind dominated hot medium be- 
tween 1.5 and 3.5 kpc, we again assume that the escap- 
ing positrons go into both the overlying halo and the 
photoionizcd, WI and HII phases in the disk. For al- 
though the upper ends of these flux tubes are blown out 
into the overlying halo, so that the positrons could all 
escape, their roots feed into the surrounding HII and 
WI envelopes that separate the HT from the HI. Thus 
within large uncertainties, we also assume the nominal 
case where escape up into the halo or down into the sur- 
rounding photoionized gas is equally probable (~ 46% 
each) , and that of those entering the photoionized phases 
the WI and HII are also equally likely (~ 23% each). 
As we show below, however, we expect that only 66% 
of the positrons that enter either the WI or the HII, 
slow down and annihilate before there, while the remain- 
der escape into the halo, where they annihilate. Thus, 
we expect that the effective net propagation factors are 
Pht-.wi ~ 0.15 ± 0.07, Pht-.hii ~ 0.15 ± 0.07, and 
Pht-.h ^ 0.70 ± 0.35, including Pht-.ht- 

5.2. Warm Neutral Medium (HI) 

Studies by Wolfire et al. (2003) find that over most 
of the Galactic disk > 3.5 kpc neutral HI gas resides in 
two stable, co-existing phases: a cold (^ 100 K) neutral 
medium and a warm (< 10^ K) neutral medium, but 
the filling factor of cold HI is negligible (Ferriere 1998). 
Further, Wolfire et al. found relative electron densities 
in the predominantly neutral HI, nni/rii ~ 1/xe ~ 50 
which greatly exceeds fcrit of equation (7), so turbulent 
Alfvenic cascades are damped by ion-neutral friction (e.g. 
Kulsrud & Pearce 1969). Thus we expect the positrons 
to stream isotropically along magnetic flux tubes at Pc/2- 

The mean HI gas density seen by the positrons escaping 
into the interstellar medium from SNIa is 0.13 cm~^, 
using the HI gas density distribution from Ferriere (1998) 
applicable at > 3.5 kpc and for our modeled SNIa and 
SNIp spatial distribution, equation (3). Thus the slowing 
down time, tgd ~ 7 x 10^ yr from equation (5), and in 
that time they are spread both directions along a flux 
tube over a slowing down length ^,5^ ~ l3c/2 ~ 90 kpc. 

We estimate the flux tube lengths in the HI in the 
Galactic disk from the mean separation between super- 
bubbles, since are the primary agents for disrupting the 
interstellar magnetic field, which otherwise tends to lie 
in the plane of the disk. From an analysis of the Galac- 
tic distribution of OB associations by McKee & Williams 
(1997), we estimate that the mean distance between su- 
perbubbles in the Galactic disk, s ~ 0.42e~'^/'^'' kpc 
as a function of galactocentric distance with Hn ^ 3.5 
kpc. Thus, the flux tube length Ib ^ s may range from 
only 1.4 kpc in the inner part of the warm disk at R 
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of 3.5 kpc to about 4 kpc at the solar radius, giving a 
mean Ig of about 3 kpc. Since these distances are far 
less that the slowing down length, we would expect that 
only a very small fraction (~ Is/'i.lsd ~ 3/180 ^ 2%) of 
the positrons, born in the HI, slow down and annihilate 
there. 

Thus, we expect that nearly all (^ 0.98J of the SNIa 
positrons, streaming along the magnetic flux tubes, escape 
from the HI into the overlying halo or the neighboring 
warm ionized HII and WI, going into the latter phases 
with equal likelihood. As we noted above, we expect that 
about 66% of the positrons entering either the WI or the 
HI, slow down and annihilate there, while the remainder 
escape into the halo and annihilate. Therefore with large 
uncertainties, we take for the net factors, Phi-.hi ~ 0.02, 
and Phiwi ~ 0.16 ± 0.08, Phihii -- 0.16 ± 0.08, and 
P///:H ~ 0.66 ± 0.33. 

5.3. Photoionized Medium in the Galactic Disk 

OB associations ionize the interstellar medium, and 
create localized, high-density (HII) regions as well as an 
extended, diffuse warm ionized medium (WIM) (e.g. Mc- 
Kee & Williams 1997). Ionizing radiation from such clus- 
tered association stars, which are embedded initially in 
their parent molecular cloud, heat as well as ionize the 
surrounding dense molecular gas. Ultimately such ion- 
ized gas expands, streams away from the molecular cloud, 
evacuating large cavities, finally destroying the molecu- 
lar cloud (Tenorio-Tagle 1979; Whitworth 1979). Thus 
OB associations create dense (> 1000 cm^^), compact 
(< 1 pc) cores of photoionized gas, observable as lumi- 
nous thermal radio sources and visible nebulae, but ^ 
2/3 of the ionizing radiation emitted by massive associa- 
tion stars escapes these dense cores and create extended 
(> 100 pc), more tenuous (< 10 cm~^) envelopes of pho- 
toionized gas (Anantharamaiah 1985; Heiles, Reach, & 
Koo 1996; McKee & Wilhams 1997). Following McKee 
& Williams we model Galactic-disk HII regions as two 
components: the dense central cores surrounding the par- 
ent OB associations and the extended HII envelopes. As 
noted by McKee & Williams the tenuous HII envelopes 
completely dominate the HII filling factor, and also about 
half the volume of an HII region is filled with a hot 
{geqKf' K) tenuous, superbubble plasma. 

Analyses of pulsar dispersion measures and diffuse, op- 
tical line emissions require the existence of a widespread 
(volume filling factor ^ 20%), warm (~ 8000 K), ion- 
ized medium (WI) of a modest density (^ 0.1 cm~'^) 
distributed throughout the Galactic disk and lower halo 
(Reynolds 1984, 1991; Tielens 2005). The extended WIM 
has a large scale height, ~ 1 kpc, transverse the Galactic 
plane (Reynolds 1991; Taylor & Cordes 1993), compared 
to ^ 0.05 kpc, the scale height of O stars (Mihalas & 
Binney 1981). However, a dilute 0-star radiation field 
seems to explain many of WI line features (Mathis 1986; 
2000). Thus it seems likely that a modest fraction, ~ 
0.15, of Lyman-continuum radiation emitted by the pla- 
nar OB associations escapes from cloud-debris HII re- 
gions to ionize the high-latitude WI (e.g. Brcgman & 
Harrington 1986; Mathis 1986; Dove & ShuU 1994; Dove, 
ShuU, & Ferrera 2000; Mathis 2000; Tielens 2005). Addi- 
tionally, Slavin, McKee, & HoUenbach (2000) found that 
the cooling of shock-heated gas, created by old supernova 
remnants located at high Galactic latitudes, creates a sig- 



nificant fraction of WI at these large kpc heights above 
the Galactic plane. 

5.3.1. Warm Ionized Medium (WI) 

Based on thermal pressure balance and the assumption 
that the WI is highly ionized, the mean electron density 
in the WI, n^ ~ 0.26 cm"^ (Kulkarni & Heiles 1987). The 
neutral fraction of the WI is not well constrained obser- 
vationally by optical-line diagnostics, but neutral frac- 
tions ~ 0.1 seem likely (Tielens 2005). However, analy- 
ses of pulsar radio signals provide a stronger constraint 
on the WI ionization fraction. Radio-wave scatterings 
increase roughly linearly with dispersion measures ^ 80 
pc cm~^, a primary WIM tracer (Cordes et al. 1991). 
Thus Cordes et al. ascertained that the thermal elec- 
trons, which for create the pulsar dispersion measures, 
also generate corresponding radio scatterings. Since such 
radio scatterings are created by electron density fluctu- 
ations on spatial scales, ~ 10^ — 10^^ cm (Cordes et al. 
1991), significantly smaller than the mean free path for 
proton-HI collisions, ^ 5 x 10^'^ cm in unit-density media, 
means that efficient turbulent radio scattering exists in 
WI and, thus, HI fraction must be less than the critical 
value, equation (7). Consequently, we assume that in WI 
Xe ^ 0.9. 

A lower limit to the typical thickness of such WI 
regions is the mean free path of 13.6 eV photons, 
(nHiani)^^, where ani ~ 3 x 10~^^ cm^ (Tielens 2005) 
and the neutral HI density, nni ^ (1 — Xe)ne/xe- Tak- 
ing a nominal value of x^ ^ 0.9 and n ~ Ue/xe ~ 0.3 
cm~^ gives a typical WI thickness, Alwi ~ 40 pc, and 
a positron slowing down time tgd ^ 1.6 x 10^ yr from 
equation (5). 

Recent simulations (Esquivel et al. 2006) of the turbu- 
lent mixing generated by shear flow instabilities in the 
WI layers between the HT and HI suggest an outer tur- 
bulent scale lo ^ 10 pc and (m^)!/^ ^ 20 km s^^. As- 
suming Bo - 4^G, {SBl^/^/Bo - 1, and Twi ^ 8000 
K, we find from equation (6) that A|| ^ 0.2 pc. Thus in 
the slowing down time of 1.6x10^ yr we would expect 
the positrons to be distributed both directions along a 
fiux tube over a distance Isd ^ (2A||/3ct,5d/3)^/^ ^ 65 pc 
in both directions along the flux tubes. 

Again assuming that the mean length l'^ of nominally 
parallel flux tubes through these thick turbulent phases 
is equal to the typical thickness of the region, Alwi ~ 40 
pc, we expect that with the turbulent scale of ^ 10 pc, 
the mean meandering tube length Ib ~ I'b/^o ~ 40^/10 ^ 
160 pc. The probability that those positrons initially 
born in the WI, slow down and annihilate there before 
they escape into the neighboring HI and HT is ~ Ib/{Ib+ 
Isd) ~ 160/225 ^ 71%. We expect that the escaping 
positrons enter into the HI and HT with equal probability 
and end up annihilating in the WI, HII or halo in the 
same relative proportions as those born in the HI and 
HT. Thus the effective propagation fractions for those 
born in the WI are Pwi-.wi ^ 0.76 ± 0.38, Pwi-.hii ^ 
0.05 ± 0.02, and Pwi-.h ~ 0.19 ± 0.09. 

On the other hand, the positrons, born in the HI and 
HT, that diffuse into the WI from are effectively reborn 
in the WI at their first scatter just inside the bound- 
ary, and since Isd < Ib, we expect that ^ 0.5 of them 
slow down an annihilate in the WI, while the other half 
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escape equally into the adjacent HI and HT. We also ex- 
pect some of the escaping positrons to subsequently enter 
other WI phases and annihilate there, or in the HII or 
halo in the same relative proportions as those born in the 
HI and HT. Thus we expect that for the positrons that 
diffuse into the WI, their effective propagation fractions 
are Pwi-.wi ~ 0.58 ± 0.29, Pwi-.hii ~ 0.08 ± 0.04, and 
Pwi-.H - 0.34 ±0.17. 

5.3.2. HII Regions 

Although the direct positron production by supcrnovae 
in the HII regions is modest because of their small filling 
factor, a large fraction of positrons made in the neighbor- 
ing HT and HI escape into it. Thus, it becomes the ma- 
jor site of positron annihilation in the disk, because of its 
high density (~ 10 times that in the HI), which reduces 
their slowing down time, so that most of the positrons 
slow down and annihilate before they can escape. 

McKee & Williams (1997) derived the Galactic volume 
filling factor of the HII regions, as a function of galacto- 
centric radius, R, assuming that these HII regions repro- 
duced the observed Galactic plane dispersion measures. 
Further, they assumed the electron density is constant 
across the Galactic disk, and all the i?-dependence re- 
sides in the filling factor, fHii{R)- They ascertained 
that fmAR) = fHii{0)e-"/"'^, where Hr = 3.5 kpc, 
and in such HII regions rig = 2.6 cm ■^, and fniii^) = 
0.13. In these HII regions, the positrons slow down in 
tsd ~ 1-7 X 10^ yr from equation (5). 

The typical thickness an HII envelope can be deter- 
mined from the average dispersion measure contributed 
by the walls of a single large superbubble HII envelope, 
ADM - 70 cmr^ pc (McKee & WiUiams 1997). For a 
mean density Ue ~ 2.6 cm^^, this corresponds to a typ- 
ical wall thickness, AIhjj ^ O.bADM/ne ~ 13 pc. So, 
again assuming nominally parallel flux tubes, the mean 
length through one of these highly turbulent walls l'^ is 
also ~ 13 pc. Although the nature of turbulent processes 
in the HII is poorly known, Haverkorn et al. (2004) have 
determined from the analysis of rotation measures that 
the outer scale of turbulent flows in disk HII regions, lo, 
is relatively small, ~ 2 pc. From this, we then estimate 
the mean meandering tube length Ib ^ I'b/^o ^ 13^/2 ^ 
85 pc. 

To compare that with the expected slowing down 
length Isd, we assume Bo ~ ifiG and Thii ~ 8000 K, 
and estimate the diffusion mean free path in the HII 
using the cooling constraint, as we have done in other 
ionized phases. Following Tielens (2005), we find that 
the cooling time is short, tcooi ~ 1-1 x 10^ yr. So tak- 
ing lo ~ 2 pc and the constraint that the dissipation 
rate of the turbulence must be less than the cooling rate 
of the plasma, we find that the root-mcan-square tur- 
bulent velocity, (u^)^/^ < 23 km s^^, which is greater 
than the Alfven speed in this plasma, Va ^ 5 km s"'^, so 



2U/2 



Bo, and from equation (6) wc thus find that 



(SB') 

A|| - 0.1 pc. 

In the slowing down time of ~^ 1.8 x 10^ yr, we would 
thus expect the positrons to be distributed over a length 
led ~ (2A||/3ctsd/3)^/^ ~ 15 pc in both directions along 
a flux tube. Thus we expect from equation (8) that 
-- h/ilB + Isd) ~ 85/(85 + 15) - 85% of the positrons 
born in one of these HII envelopes, would slow down and 



annihilate there before they escape into the neighboring 
HI and HT. We again assume that half of the escap- 
ing positrons go into each of these two phases, and that 
from there some diffuse into other HII and WI regions 
or into the halo and annihilate in the same proportion 
as those initially born in the HI and HT. Thus the ef- 
fective propagation fractions for those born in the HII 
are Phiihii -- 0.87 ± 0.43, Phiiwi - 0.03 ± 0.01, and 
i^H//:H-^ 0.10 ±0.09. 

Again, however, most of the positrons in the HII were 
born in the HI and HT, and were reborn in the HII within 
a scattering mean free path of the boundary. As we have 
discussed above, with Isd < Ib, we expect that ^ 0.5 of 
the positrons slow down an annihilate in the HII, and 
the other half escape equally into the adjacent HI and 
HT. From there, we again expect them to eventually an- 
nihilate in the HII, WI or halo in the same relative pro- 
portions as those born in the HI and HT. Thus, as with 
the WI discussed above, we expect that for the positrons 
that diffuse into the HII, their effective propagation frac- 
tions are Phii-.hh -- 0.58 ±0.29, Phii-.wi -^ 0.08 ±0.04, 
and Phii-.h ~ 0.34 ±0.17. 

5.4. Disk Annihilation Flux Asymmetry 

The asymmetry in the 511 keV line flux (Weidenspoint- 
ner et al. 2008b) from the inner disk provides further evi- 
dence of positron propagation and annihilation in the HII 
and WI. We show that this flux asymmetry can be fully 
accounted for by the apparent asymmetry of positron an- 
nihilation in the inner spiral arms, as viewed from our So- 
lar perspective, without any asymmetry in the positron 
production. 

Modeling the inner Galactic disk in two compo- 
nents, Weidenspointner et al. (2008b) find a best-fit 
SPI/INTEGRAL 511 keV line flux of (4.3 ± 0.5) x lO"'^ 
photons cm~^ s"'^ from the component between —50° < 
/ < 0°, compared to a flux of (2.4 ± 0.5) x lO""* pho- 
tons cm^^ s^^ from that between +50° > Z > 0°. This 
amounts to an excess of 80 ± 40% in the 511 keV flux 
from negative longitudes compared to that from positive 
longitudes, and it amounts to an asymmetry of about 
20±8% in the total observed disk flux (Weidenspointner 
et al. 2008a). 

The asymmetry in the Galactic spiral arms, as viewed 
from our solar perspective, can be seen in Figure 4, show- 
ing the current mapping of the four spiral arms deter- 
mined from the study by Vallee (2005a, b). As summa- 
rized in Table 4, we expect that the positron annihilation 
in the disk occurs almost entirely in the warm ionized HII 
and WI phases, which are concentrated along the spiral 
arms (e.g. Drimmel & Spergel 2001). Thus, we assume 
that equal fractions of the positrons produced at any par- 
ticular radius propagate to and annihilate in the ionized 
gas in their nearest spiral. 

We estimate the positron annihilation rate per unit 
length along each of the four arms, using the radial dis- 
tributions of the positron production surface density in 
the disk as a function of Galacto-centric radius, R, given 
in equation (3), multiplied by T:RdR/2 and by the disk 
annihilation/production fraction (Table 4). We further 
weight the annihilation distribution along the flux tubes 
by a unit normalized surface density distribution of in- 
terstellar gas, since positrons propagating along the flux 
tubes are more likely to encounter the dense ionized gas. 
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where they annihilate, in denser regions. We also in- 
clude the positrons escaping into the arms from the outer 
bulge. In accordance with observations (e.g. Drimmel & 
Spergel 2001) that the Sagittarius arm is only half as 
dense and active as the others, we also reduce it accord- 
ingly. 

At every point on each arm, we then multiply the dif- 
ferential annihilation rate by the 75ii/e+ ratio to deter- 
mine the differential 511 keV line luminosity and then 
divide by 4:TrdQ, where d is the distance from the sun, 
to calculate the expected differential 511 keV line flux 
at the earth. Finally, integrating along each arm out to 
±50° Galactic longitude within ±10° latitude to compare 
with the SPI/INTEGRAL model, we find asymmetric 
511 keV line disk fluxes of (3.8±1.0) x lO""* photons cm"^ 
s^^ from the negative longitudes, and (2.3 ± 0.6) x 10^^ 
photons cm~^ s~^ from the positive, and a negative-to- 
positive longitude flux ratio of 1.6±0.4. These values 
are all quite consistent (Table 7) with the corresponding 
best-fit (Weidenspointner et al. 2008b) SPI/INTEGRAL 
line fluxes of (4.3 ± 0.5) x 10""' photons cm~^ s~^ and 
(2.4 ± 0.5) X 10~^ photons cm~^ s~^, and a negative-to- 
positive flux ratio of 1.8±0.4. 

Thus we show that the observed (Weidenspointner et 
al. 2008b) SPI/INTEGRAL asymmetry m the 511 keV 
line flux from the inner disk can be fully explained by the 
apparent asymmetry in the inner spiral arms as viewed 
from our Solar perspective, which leads to different line- 
of-sight intensities, without any difference in the positron 
annihilation or production rates. This asymmetry in the 
annihilation radiation from the Galactic disk is a product 
of the Norma arm dominating the negative longitudes due 
to its proximity to the inner region of the Galaxy, and 
the weakness of the Sagittarius- Garina arm due to its 
greater distance from the Galactic center from which the 
annihilation rate falls off exponentially from the Galactic 
center and its lower density /activity (Drimmel & Spergel 
2001) compared to the other 3 arms. 

5.5. Disk Positron Annihilation >3.5 kpc 

We now estimate the total outer disk annihilation rate 
from the contributions of the various phases and the 
propagation fraction within and between them. These 
are summarized in Table 4. 

As we have discussed above, we see that more than 
half of the positrons, produced in the dominant phases 
of the interstellar medium in the disk, the HI and HT, 
escape into adjacent phases before they can slow down 
and annihilate. Only those produced in the HII and WI 
phases have a large chance of annihilating before they 
escape, and the rest either end up in the HII, the WI, or 
the halo. This result is quite robust, even allowing the 
large uncertainties in the relative likelihood of positrons 
following magnetic flux tubes that lead into these three 
annihilation traps. 

Thus, we expect that the total annihilation rate in the 
disk beyond 1.5 kpc is Ad ~ (0.59 ± 0.11) x lO''^ e+ 
s^^. As we discuss in detail in the summary (§ 7), this is 
quite consistent with the best-fit disk annihilation rate 
of inferred from the SPI/INTEGRAL analyses (Weiden- 
spointner et al. 2007, 2008a). 

6. POSITRON PROPAGATION & ANNIHILATION IN THE 
GALACTIC HALO 



Lastly we consider the halo, which for comparison with 
the SPI/INTEGRAL analyses we take to lie beyond 1.5 
kpc from the Galactic center and more than 0.4 kpc 
above the disk. Although it has no direct positron pro- 
duction by supernovae, a large fraction of positrons made 
in all the underlying phases of the bulge and disk can es- 
cape into it. Thus, it too becomes a major site of positron 
annihilation, because its pervasive hot tenuous plasma, 
coming primarily out of the outer bulge and disk super- 
bubbles (e.g. Norman & Ikeuchi 1989), may be turbu- 
lent enough to produce a scattering mean free path much 
smaller than its size, so that most of the positrons that 
enter it slow down and annihilate before they can escape 
on out of the Galaxy. 

From the positron production in all of the regions and 
phases and their propagation fractions, discussed above, 
we estimate the positron escape into the halo and their 
annihilation there to be, 



Ah ^ PH:H{QBiPBi:H + QsmPBrn-.H 

-QboPbo-.h + QhtPht-.h + QhiPhi-.h 

+ QwiPwi:H + QhIiPhII:h), 



(17) 



Extensive observations of X-ray emission from the 
Galactic halo (e.g. Snowden et al. 1998; Pietz et al. 
1998) have shown that hot tenuous plasma filling it is 
essentially indistinguishable from that in the outer bulge 
and superbubbles from which it is derived. Thus we 
assume a similar T of 1x10® K, a density of 2x10""^ 
cm^^, and magnetic field, _Bo ^ 4 jiG. This plasma 
has a scale height, zxr, ~ 4.4 kpc, and a radial scale 
length of 15 kpc. As discussed above for comparison 
with SPI/INTEGRAL analyses, we define the halo as 
that region beyond a radius 1.5 kpc and a z > 0.4 kpc, 
including within it the outer bulge. 

The nature of turbulence in the halo is very uncertain, 
but we assume that a significant level of MHD turbu- 
lence can be maintained by erupting superbubbles. Using 
the parameters for that hot plasma, an outer turbulence 

scale, lo ~ 75 pc and {B'^Y^'^/Bo < 0.5, we also find the 
expected scattering mean free path parallel to the mean 
magnetic field, Aii ^ 14 pc from equation (6). 

Assuming a minimum flux tube escape length Ig ^ 
10 kpc, roughly 2 times the scale height, then the ex- 
pected positron escape time tesc ^ 200 Myr. By com- 
parison the positron slowing down time, tsd ~ 23 Myr 
from equation(5), and in that time positrons travel Igd ^ 
(2/3ctsdA|i/3)^/^ ^ 7 kpc. Since all of the positrons enter 
the halo from below, rather than being born throughout, 
wc expect that they will all slow down and annihilate 
well before they can escape, since Isd < Ib- Thus, for all 
of the positrons that enter the halo, we expect Ph-.h ~ 1- 
Because the properties of the halo are assumed to be the 
same and those of the outer bulge and superbubble HT, 
which are the major sources of positrons that diffuse into 
it, there is no effective boundary between them and no 
"reflection" of those positrons from the halo. Moreover, 
as we discussed above, all positron annihilation in the 
outer bulge and superbubble HT is also included in the 
halo, because of their large scale sizes. 

From the estimates of the positrons entering the halo 
from the various underlying bulge and disk phases, sum- 
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marized in Tables 4 and 5, we expect a total positron 
annihilation rate in the halo beyond 1.5 kpc of ~ (0.82 ± 
0.23) X 10''3 e+ s^i. As we show in Table 6 and dis- 
cuss further in the summary, this annihilation rate gives 
an average 511 keV line flux from the halo of ^ (1.58 ± 
0.37) X 10~^ cm"^ s~^, depending on how much refrac- 
tory dust survives in the halo. This flux together with 
the bulge flux of ~ (0.74 ± 0.14) x 10"^ cm"^ g-i gj^gg 
a combined bulge and halo flux of ~ (2.32 ± 0.40) x 10~^ 
cm~^ s~^ (see Table 7), which is quite consistent with 
that of ^ (2.14 ± 0.12) x 10~^ cm'^ s"\ determined 
from the SPI/INTEGRAL best-fit flux model (Weidens- 
pointner et al. 2008a). 

As we discuss further in the summary, not only the 
511 keV line intensity but also the Galactic averages of 
both its broad/narrow line flux ratio and the positron- 
ium fraction are very sensitive to the mean dust content 
of the hot halo plasma. For as Jean et al. (2006) have 
shown and we list in Table 6, if the refractory elements 
are all in dust in the hot plasma, only ^ 18% of the re- 
sulting positron annihilation occurs via positronium for- 
mation, and fully ^ 90% of the 511 keV line emission is 
narrow (FWHM ^ 2 keV) while the remaining 10% is 
very broad (FWHM ^11 keV). This gives a halo flux of 
~ (1.60 ±0.37) X 10^'^ cm^^ s~^ in the narrow line alone 
which in turn gives a combined bulge and halo flux of 
~ (2.34 ± 0.40) X 10"-'' cm"2 s"\ which is likewise quite 
consistent with that from the SPI/INTEGRAL best fit 
above. Combined with disk flux, that would give a total 
Galactic bulge, disk and halo flux of - (3.25±0.42)xl0"3 
cm"^ s~^ which is also quite consistent with the com- 
bined best-fit flux of - (2.87 ± 0.25) x lO'^ cm'^ s'^, 
from the SPI/INTEGRAL analyses (Weidenspointner 
ct al. 2008a). This also gives Galactic average val- 
ues of the positronium fraction of only ~ 0.66 ± 0.12 
and the broad/narrow 511 keV line flux ratio of only 
~ 0.17 ± 0.04, for which the best-fit SPI/INTEGRAL 
values have not yet been published. 

If on the other hand, the grains have all disintegrated, 
the halo flux would still be also as large, ^ (1.41±0.31) x 
10~^ cm~^ s~^, but it would all be in the very wide (^ 
11 keV) 511 keV line, so the Galactic average value of 
the broad/narrow 511 keV line flux ratio would jump to 
~ 1.30 ± 0.29, and the positronium fraction would in- 
crease only slightly to ^ 0.74 ± 0.13. The actual dust 
content of the halo undoubtedly lies in between and the 
SPI/INTEGRAL observations can provide a unique mea- 
sure. 

This is also consistent with that implied by the 
SPI/INTEGRAL best-fit flux model of Ah ~ (0.65 ± 
0.11) X 10''^ e+ s-i (Weidenspointner et al. 2008a), scaled 
to a galacto-centric distance of 8 kpc and (e+/75ii) ~ 
0.65±0.07 expected (Jean et al. 2006) in the halo plasma, 
as discussed above. 



7. SUMMARY 

The recent analyses by Weidenspointner et al. (2007, 
2008a) of the SPI/INTEGRAL measurements of the 
Galactic 511 keV positron annihilation radiation lumi- 
nosity suggest that the best flt to the observations con- 
sists of three basic components: 1) a spherical bulge of ~ 
1.5 kpc radius, 2) a thick (± 400 pc) Galactic disk, and 3) 
a spherical halo extending beyond 1.5 kpc. Their most 



recent best-fit model (Weidenspointner et al. 2008a), 
scaled to a galacto-centric distance of 8 kpc, gives 511 
keV line luminosities of is ^ (0.57±0.08) x 10''^ photons 
s-i from the bulge, Ld ^ (0.40±0.06) x lO''^ photons s'^ 
from the disk, and a net Lh ^ (1.0 ±0.1) x 10"''' photons 
s^^ from the halo beyond the bulge, for a total Galactic 
511 keV luminosity L - (2.0 ± 0.2) x lO'*^ photons s~^ 

This gives a positron annihilation luminosity 
bulge/disk ratio of ~ 1.4±0.3. But since this is a 
factor of 4 higher than the expected positron production 
ratio, inferred from the stellar bulge/disk ratio of the 
assumed Galactic supernova sources, various exotic 
new sources have been invoked to explain the large 
discrepancy. 

Here we have shown, however, that the measured 511 
keV luminosity ratio can be fully explained by positrons 
from the decay of radionuclei made by explosive nucle- 
osynthesis in supernovae, if the propagation of these rel- 
ativistic positrons in the various phases of the interstellar 
medium is taken into account, since these positrons must 
first slow down to energies < 10 eV before they can anni- 
hilate, and if the geometry of the sources is also consid- 
ered. Moreover, as Jean et al. (2006) have shown, with- 
out propagation none of the proposed positron sources, 
new or old, can explain the two fundamental properties 
of the Galactic annihilation radiation: the fraction of the 
annihilation that occurs through positronium formation 
and the ratio of the broad/narrow components of the 511 
keV line. 

7.1. Positron Production 

First, we have shown, using recent estimates (Cap- 
pellaro, Evans & Turatto 1999) of the SNIa and SNIp 
rate in our Galaxy of ~ 0.40±0.16 SN per 100 years, 
that the positrons resulting from the /3+-decay chains of 
the radioactive nuclei, ^^Ni, ''^Ti, and ^^Al, produced in 
these supernovae can fully account for the total Galactic 
positron production rate of ^ (2.4 ± 0.6) x 10^'^ e+ s~^, 
imphcd by the best-fit SPI/INTEGRAL analyses (Wei- 
denspointner et al. 2008a). 

The total Galactic production rate of positrons from 
long lived (1.04 Myr meanhfe) decay of ^®A1 of ~ (0.3 ± 
0.1) X 10'''^ e+ s~^, produced by massive stars, is deter- 
mined directly from the SPI/INTEGRAL measurements 
(Dichl et al. 2006) of the Galactic luminosity of the 1809 
keV line emission which accompanies the decay. The 
Galactic positron production rate the ^'^Ti decay chain 
(89 yr meanlife) of - (0.5 ±0.2) x lO^^ g-H g-i ^an be de- 
termined from the relative Galactic abundances (Lodders 
2003) of ^^Ca and ^^Fe, which are primarily produced by 
'''^Ti and ^^Ni decays, using the calculated SNIa rate and 
^®Ni yield, assuming that these supernovae presently pro- 
duce about half of the ^®Fe and core collapse supernovae 
produce the other half. The Galactic positron produc- 
tion from the much shorter lived (111.4 day meanlife) 
decay chain ^^Ni, however, depends on an additional fac- 
tor, since most of these decay positrons slow down and 
annihilate in the dense ejecta and only a fraction, /sg 
survive into the interstellar medium. We showed that 
all of the additional positron production, amounting to 
~ 1.6 X 10^^ e+ s^^ , can be accounted for by the ^^Ni de- 
cay chain with a positron survival fraction in SNIa ejecta 
of /56 ^ 5±2%. This is the time-integrated mean of ^ 
5% calculated (Chan & Lingenfelter 1993 Fig. 3) from 
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the standard deflagration models of SNIa (Nonioto ct al. 
1984, W7), assuming combed out magnetic fields, and it 
is also consistent with the mean fraction of ~ 3.5 ± 2% 
inferred (Milne, The & Leising 1999) from SNIa hght 
curves at late times. 

Although the uncertainties in both the SNIa rate and 
the normalized survival fraction are large, ^ 40 to 50 %, 
they are necessarily anticorrelated. Therefore their prod- 
uct yields an uncertainty of only ^ 25% in the production 
rates. 

Using these positron production rates and the Galac- 
tic bulge and disk radial distributions of the SNIa, SNIp, 
and massive stars, gave an expected positron production 
of (0.31 ± 0.07) X lO'*^ e+ s'^ in the inner bulge < 0.5 
kpc, (0.37 ±0.09) X 10"*^ e+ s^^ in the middle bulge from 
0.5 to 1.5 kpc, (0.56 ± 0.14) x lO'^^ e+ s^^ in the outer 
bulge from 1.5 to 3.5 kpc, and (1.16±0.22) x lO'*^ e+ s'^ 
in the disk beyond 3.5 kpc. For a positron annihilation 
bulge within 1.5 kpc used in SPI/INTEGRAL observa- 
tion analyses, that gave a positron production bulge/disk 
ratio of 0.4 ±0.1 within 1.5 kpc. But, as we have shown, 
the expected annihilation ratio is very different, because 
the interstellar medium and the resulting positron prop- 
agation in these regions is also very different. 

7.2. Positron Propagation 

We have modeled the properties of the various regions 
and phases of the interstellar medium to determine the 
propagation of the positrons as they slowed down and an- 
nihilated. In ionized plasmas we assumed that diffusive 
propagation was controlled by resonant scattering of the 
relativistic positrons by MHD waves at their cyclotron 
radius. These waves are generated by the cascade of 
the larger scale magnetic turbulence down to sufficiently 
small scales which scatter the positrons and they diffuse 
along magnetic flux tubes. Such propagation has been 
observed in interplanetary plasmas with MeV electrons 
and we use a phenomenological model based on these 
studies. In neutral, or at least nearly neutral, phases of 
the interstellar medium, however, this cascade is damped 
by ion-neutral collisions, and the positrons are expected 
to stream along the magnetic flux tubes with an isotropic 
pitch angle distribution. 

From these propagation processes, we have calcu- 
lated the propagation, slowing down and annihilation of 
positrons formed in the various regions and phases of the 
interstellar medium described in Tables 1 and 3. From 
the properties of the medium in each phase, we deter- 
mined the propagation mode, and calculated the diffu- 
sion mean free path in the undamped cascade mode, or 
the streaming velocity in the damped mode. From an 
estimate of the mean length of the magnetic flux tubes 
within a phase, compared to the tube length over which 
the positrons would propagate before they slowed down, 
we estimated the fraction of positrons expected to slow 
down and annihilate in that phase before they escape. 
We then estimated the relative fractions of the escaping 
positrons that go into each of the adjacent phases by sim- 
ple geometric arguments. We have assumed large uncer- 
tainties of ± 50% in the estimated propagation fractions 
between phases, which produce comparably large uncer- 
tainties in the annihilation rates in each phase, although 
their sums within the bulge and disk are more tightly 
constrained. These propagation fractions between var- 



ious components and the resulting distribution of their 
annihilation rates are presented in Tables 2 and 4. 

7.3. Positron Annihilation Rate & Flux 

The expected positron production and annihilation 
rates in the different regions and phases of the interstel- 
lar medium are summarized in Tables 5 and 6. Finally, 
the expected positron annihilation rates and 511 keV line 
properties are compared with SPI/INTEGRAL analyses 
in Table 7. 

We show that ~ 52±7% of the all the Galactic 
positrons are expected to be produced within the bulge < 
3.5 kpc, of which ^ 28±6% are born within the positron 
bulge < 1.5 kpc. From the relative filling factors of 
the various gas and plasma phases in the bulge, we ex- 
pect that essentially all of these positrons are born in 
the hot fully ionized plasma. From their propagation 
fractions we also expect that ~ 80% of them diffuse 
into the warm HII and cold HI envelopes of molecular 
clouds that lie within 1.5 kpc, where they slow down 
and annihilate, while the remaining ~ 20% escape into 
the halo and disk beyond. This propagation thus re- 
sults in a positron bulge annihilation rate <1.5 kpc of 
Ab ~ (1.00 ± 0.20) X 10^3 g-H g-i^ 

Of the positrons that are produced in the Galactic disk 
beyond 3.5 kpc, we expect from the relative filling factors 
that that about 70% of the disk positrons are born either 
in the ubiquitous warm neutral HI medium of the disk, 
or in the hot tenuous plasmas of superbubbles, in both 
of which the conditions are such that positrons rapidly 
stream or diffuse along magnetic flux tubes. Some seg- 
ments of these flux tubes thread through the warm ion- 
ized gas and HII shells that separate the warm neutral 
gas from the hot tenuous plasmas, while other segments 
have been blown up into the halo by the hot expand- 
ing superbubbles (e.g. Parker 1979), and, since there is 
no preferred direction, the positrons diffusing or stream- 
ing along them should have roughly equal likelihood of 
either going into the warm denser ionized shells, where 
they mostly slow down and annihilate, or out into the 
halo. From the estimated propagation fractions, assum- 
ing large uncertainties, we flnd that nearly half of the 
positrons enter the ionized shells and annihilate in the 
disk, and that other half escape into the overlying halo, 
giving a total annihilation rate in the disk beyond 3.5 
kpc is Ad ~ (0.59 ± 0.11) x lO*^ g+ g-i 

So even though roughly equal numbers of positrons are 
born in the bulge within 3.5 kpc and in the disk beyond, 
their annihilation rate is higher in the bulge than in the 
disk, because ^ 80% of those born in the interstellar bulge 
slow down and annihilate in the cloud envelopes within 
1.5 kpc, while only ~ 51% of those born in the disk slow 
down and annihilate before they escape into the halo. 

We now compare the spatial distribution of the Galac- 
tic positron annihilation radiation, expected from the 
/J^-decay chains of the radioactive nuclei, ^^Ni, ^^Ti, 
and ^^Al, produced in supcrnovae, with that observed 
by SPI/INTEGRAL. 

The most direct comparison is through the 511 keV line 
fluxes. From more than four years of SPI/INTEGRAL 
observations Wcidcnspointner et al. (2008a, model BD) 
flt the 511 kcV fluxes and luminosities to a combined 
bulge component within a spheroidal Gaussian (FWHM 
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of 1.5 kpc) distribution and thick disk component. They 
found best-fit fluxes of - (0.75 ± 0.09) x lO^-'^ photons 
cni-2 s-i from the bulge and - (0.94 ± 0.16) x 10"-'' 
photons cm~^ s~^ from the disk. They also found a best- 
fit combined halo and bulge flux of ^ (2.14±0.11) x 10~^ 
photons cm~^ s~^. 

For comparison we calculate in Table 5 the expected 
fluxes from the annihilation rates in the different re- 
gions and phases listed in Tables 2 and 4, using the 
511 keV photons yields per positron expected in those 
phases (Guessoum, Jean & Gillard 2005; Jean et al. 
2006). Thus, we expect 511 kcV line fluxes of ^ 
(0.72 ± 0.13) X 10-3 photons cm^^ g-i ^^^^ ^^^ bulge, 
~ (1.08 ± 0.16) X 10-3 photons cm^^ g-i ^^^^ ^^le disk, 
and ^ (2.37 ± 0.41) x 10-^ photons cm"^ s~^ from the 
combined halo and bulge. These expected fluxes, as sum- 
marized in Table 7, are all in excellent agreement with 
the best-fit fluxes observed by SPI/INTEGRAL (Wei- 
denspointner et al. 2008a). 

We also calculated the expected 511 keV fluxes from 
the inner bulge within 0.5 kpc and the middle bulge 
between 0.5 and 1.5 kpc, - (0.39 ± 0.10) x 10"^ ^nd 
~ (0.33 ± 0.07) X 10-3 photons cm-^ s"\ respec- 
tively. As mentioned above, these boundaries were orig- 
inally chosen for our calculations because these were the 
two regions modeled in the previous analysis by Wei- 
denspointner et al. (2007) of the first two years of 
SPI/INTEGRAL, where they found best-fit fluxes of 
- (0.38 ± 0.03) X 10-3 and ~ (0.41 ± 0.06) x 10-3 pho- 
tons cm^^ s^^ from the same regions. These expected 
and observed values are also consistent within the I — a 
uncertainties. 

We have compared with the earlier analysis, because 
the fluxes in the later Gaussian component analysis by 
(Weidenspointner et al. 2008a) overlapped in the inner 
volume which was also smaller than 0.5 kpc. The best- 
flt total fluxes from the analyses of (Weidenspointner et 
al. 2007) and (Weidenspointner et al. 2008a) within 1.5 
kpc are essentially the same, ^ (0.79 ± 0.07) x 10-3 and 
~ (0.75 ± 0.09) X 10-3 photons cm^^ s^^, respectively. 

We also compare the expected positron annihila- 
tion rates in the bulge, disk and halo, and the cor- 
responding bulge/disk ratio, with those determined 
by Weidenspointner et al. (2008a) from their bcst- 
flt SPI/INTEGRAL observations. Since they used a 
galacto-centric distance of 8.5 kpc rather than the 8 kpc, 
which we assumed, we scale our annihilation rates by 
(8.5/8)^ for comparison. In addition, since ^ 12% of 
their disk annihilation component lies within 1.5 kpc, 
we further modify our expected annihilation rates to 
compare with their model components. Thus, to com- 
pare with the Weidenspointner et al. (2008a) values, 
we calculate an equivalent SPI disk annihilation com- 
ponent spiAd ^ ^1,(8.5/8)7(1 - 0.12) - (0.76 ± 
0.13) X 10^3 e+ s-^. This is in very good agreement 
with the SPI/INTEGRAL value of (0.81 ± 0.14) x 10'*3 
e^ s~^. Similarly, we calculate an equivalent SPI bulge 
annihilation component spiAb ^ {As — A/, [0.12/(1 — 
0.12)])(8. 5/8)2 ^ (1.04 ± 0.21) x 10'^3 ^+ g-i^ ^i^^^,!^ ^g 
also in very good agreement with the SPI/INTEGRAL 
value of (1.15 ± 0.16) x 10''3 e+ s'^. 

Thus, we expect an equivalent bulge/disk ratio of the 
annihilation rate spiAb/spiAd ^ (1.04± 0.21)/(0.76± 
0.13) ~ 1.37 ±0.37, which is also in excellent agreement 



with the measured value of (1.15 ± 0.16)/(0.81 ± 0.14) - 
1.42 ± 0.32 from the SPI/INTEGRAL data analyses of 
the best-fit bulge and disk components (Weidenspointner 
et al. 2008a). By comparison we also note, however, 
that the ratio of the total bulge annihilation with 1.5 
kpc compared to that in the disk beyond is Ab/Au ^ 
(1.00±0.19)/(0.59±0.11)~ 1.69±0.4'5. This is nonethe- 
less consistent with the SPI/INTEGRAL best-fit model 
when the 12% disk component within 1.5 kpc is sub- 
tracted from their disk component and added to that of 
the bulge, - (0.94 + 0.08)/(0.66 - 0.08) - 1.79 ± 0.38. 

We also showed that positron propagation and anni- 
hilation in the HII and WI gas along the spiral arms 
account for the asymmetry in the SPI/INTEGRAL 511 
keV line flux (Weidenspointner et al. 2008b) from mod- 
eled inner disk components within Galactic longitudes 
±50° that show a best-fit excess of 80±40% in the 511 
keV line flux from the negative longitudes compared to 
that from the positive longitudes. As we have shown (§ 
4.5 and Table 7), this flux asymmetry can be fully ex- 
plained by the apparent asymmetry in the inner spiral 
arms as viewed from our Solar perspective, which leads 
to different linc-of-sight intensities, without any actual 
asymmetry in the positron annihilation or production 
rates. 

From the estimated fractions of the positrons entering 
the halo from the underlying bulge and disk, summarized 
in Table 4, we also compare the total expected positron 
annihilation rate in the halo with the SPI/INTEGRAL 
best-fit model of a combined bulge and halo component 
by Weidenspointner et al. (2008a). To do so we combine 
the expected bulge (< 1.5 kpc) and halo (> 1.5 kpc) an- 
nihilation rates from Table 5, scaled to a galacto-centric 
distance of 8.5 kpc, to give an equivalent SPI annihi- 
lation component, spiAh+b ~ {Ab ± Ah){^.^/^Y ^ 
(1.95 ± 0.28) X 10'*3 e+ s-^ From their combined halo 
and bulge fiux, discussed above, Weidenspointner et al. 
(2008a) estimate a significantly larger combined annihi- 
lation rate of - (3.13 ± 0.2) x 'l0''3 e+ s^^. 

But they use a uniform e+/75ii ^ 1.82 to convert from 
observed the 511 keV photons luminosity to an assumed 
positron annihilation rate. This conversion factor, how- 
ever, was based on a positronium fraction of ^ 97%, de- 
termined by Jean et al. (2006) for warm gas in the inner 
Galaxy. Although this would be applicable to the bulge 
contribution, it is not applicable to the hot plasma of the 
halo contribution, for which Jean et al. (2006) estimate a 
positronium fraction of only about 18% to 42%, depend- 
ing on the refractory grain abundance, which corresponds 
to a e"'"/75ii ~ 0.65 ± 0.07, as discussed above. Thus 
we use this much more appropriate conversion factor for 
flux contribution from the halo beyond the bulge, which 
from the best-flt fluxes amounts to roughly 65% of the 
combined halo and bulge flux. We estimate an annihila- 
tion rate of spiAh ^ (0.65/1.82)(3.13 - 1.15) x 10^3 _ 
(0.71 ± 0.26) X 10''3 g+ g-i jus^ foj. the halo beyond 
the bulge, which, adding back the bulge contribution, 
gives an adjusted best-fit SPI/INTEGRAL model of 
Ah+b ^ (1.86±0.26) X 10''3 e+ s^^. Those values would 
be in good agreement with the expected equivalent SPI 
annihilation components, spiAr ^ (0.91 ± 0.28) x 10^3 
and SPI Ah+b ~ (1-95 ± 0.28) x 10^3 ^+ g-i_ 

7.4. Positronium Formation 
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Significant propagation of the positrons also accounts 
for the other fundamental features of the Galactic an- 
nihilation radiation: 1) the high, observed fraction 
of positrons annihilating via positroniuni formation of 
94±7%, a weighted average of 94±6% (Churazov et al. 
2005), 95±3% (Jean et al. 2006) and 92±9% (Weidens- 
pointner et al. 2006), and 2) the observed ratio of broad 
(^ 5.8 keV) to narrow (^ 1.2 keV) 511 keV line emis- 
sion from the bulge of ~ 0.5±0.2 (Chuzarov et al. 2005; 
Jean et al. 2006), which gives a measure of the frac- 
tion of positronium that is formed by charge exchange 
with neutral HI in flight and results in broad line emis- 
sion. For even though roughly 90% of the positrons in 
the bulge (< 1.5 kpc) are either born in, or diffuse into, 
the hot plasma that fills most of the bulge, hot plasma 
can not be the site of most of the annihilation because 
neither the observed positronium fraction, nor the broad 
to narrow line ratio, are at all consistent with annihila- 
tion in that phase (Guessoum, Ramaty & Lingenfelter 
1991; Guessoum, Jean & GiUard 2005).' 

We show that the expected positron propagation, slow- 
ing down and annihilation in the bulge and disk quite 
naturally accounts for both of these observations. In par- 
ticular, we expect (Table 5) that (0.58 ±0.11) x lO''^ g+ 
s~^ annihilate in the warm ionized and HII phases of the 
disk, (0.46 ±0.08) x 10^"^ e+ s^^ annihilate in the warm 
HII phase of the bulge, and (0.48 ± 0.17) x lO'*^ e+ s^^ 
annihilate in the cold neutral phase of the bulge. In these 
phases positronium is formed ~ 89% and ~ 95% of the 
time, respectively (Guessoum, Jean & GiUard 2005; Chu- 
razov et al. 2005; Jean et al. 2006), so we would expect 
(as outlined in Table 6) a combined positronium annihi- 
lation rate of (1.40 ± 0.74) x lO''^ e+ s^^ out of a total 
positron annihilation rate of (1.53 ± 0.81) x 10^"^ e"*" s~^ 
for a total Galactic positronium fraction fpssim92±2%. 

We also expect (Table 6) that of the positronium 
formed in the warm ionized and HII phases ~ 56% is 
in the disk and ~ 44% is in the bulge, while only ^ 2% 
of that formed in the neutral phases is in the disk and ^ 
98% is in the bulge. This gives an expected positronium 
fraction of ~ 93±2% in the bulge and ~ 90±2% in the 
disk (Table 6). The mean positronium fraction of fps ~ 
93±2% of the Galactic positron bulge within 1.5 kpc is in 
very good agreement with the mean of 9^±6% measured 
by SPI/INTEGRAL (Churazov et al. 2005; Jean et al. 
2006; Weidenspointner et al. 2006) from the inner ±20° 
of the Galaxy. 

7.5. Broad/Narrow 511 keV Line Ratio 

We also consider the ratio of the broad and narrow 
components of the 511 kev line fluxes, since the bulge 
and disk values are expected to be quite different (Table 
6), and compare the expected values with those measured 
by SPI/INTEGRAL (Table 7). 

In the warm ionized and HII phases, essentially all 
(~99%) of the positronium is formed by radiative combi- 
nation after the positrons are thermalizcd and its annihi- 
lation produces narrow (^1.2 keV FWHM), 27 511 keV 
emission 25% of the time and 37 continuum emission the 
rest of the time. In the cold neutral phases essentially all 
of the positronium is formed by charge exchange in flight 
producing the broad (^ 5.8 keV) doppler shifted 511 kcV 
line, while in the warm neutral phase, which accounts for 
a negligible ~ 2% of the neutral phase positronium for- 



mation, the fraction of in-flight formation is ^ 70% for 
the typical weakly ionized {xe ^ 2%) gas considered here. 

Thus, we would expect a total Galactic positron an- 
nihilation rate of (0.47 ± 0.12) x lO'^^ e+ s'^ from 
positronium formed in-flight and an annihilation rate of 
(0.93±0.11) X 10^'^ c+ s^^ from that formed thermally by 
radiative combination. With positronium producing two 
511 keV photons in 25% of its annihilations, this gives 
a broad 511 keV Hue luminosity of (0.24 ± 0.06) x lO''^ 
73^^ from in-flight annihilation and a narrow line lumi- 
nosity of (0.46 ± 0.06) X 10**^ 7 s"'^ from thermal annihi- 
lation. However, the residual direct thermal annihilation 
of (0.12 ±0.06) X 10**^ e+ s~^ from non-positronium pro- 
cesses also produces 2 narrow 511 keV photons from ev- 
ery annihilation with a luminosity of (0.26 ±0.12) x 10''^ 
7 s~^, which gives a total narrow 511 keV line luminosity 
of luminosity of (0.72 ± 0.14) x lO'^^ ^ g-i 

We expect (Tables 6 & 7) that of the bulge and disk 
annihilation that 68% occurs in the warm ionized gas 
versus 32% in the neutral gas. Nearly all (~ 98%) of the 
annihilation in the disk is in the warm ionized gas, while 
in the bulge we expect that 49±12% is in the ionized gas 
versus 51±14% in the neutral. The latter values are quite 
consistent with the SPI/INTEGRAL measurements from 
the bulge region, where Jean et al. (2006) find 51±3% 
in the ionized gas versus 49±3% in the neutral. 

Wc further expect (Table 6) that of the narrow line 
positronium formed in the warm ionized and HII phases 
~ 56% is in the disk and ~ 44% is in the bulge, while 
again only ~ 2% of that formed in the neutral phases is 
in the disk and ~ 98% is in the bulge. In the bulge this 
yields a broad 511 keV line luminosity of (0.23 ± 0.08) x 
10^'^ 7 s~^ from in-flight annihilation and a narrow line 
luminosity of (0.34 ± 0.06) x lO**^ 7 s~"^ from thermal 
annihilation, for a broad/narrow line luminosity ratio of 
0.68±0.26. In the disk the broad line luminosity is only 
^ 0.01 X 10*^ 7 s~^ and the narrow line luminosity is 
(0.38 ± 0.06) X 10^"^ 7 s^^, yielding a broad/narrow line 
luminosity ratio of only ^ 0.03, and the Galactic mean 
broad/narrow line luminosity ratio from the combined 
bulge and disk is 0.33±12. 

However, the apparent ratio that would be seen from 
the Earth if the full disk were viewed is quite different, 
because of different relative weighting, since the observed 
fluxes from the bulge are equal to the bulge luminosi- 
ties divided by AttR\ ~ 7.65 x 10**^ cm^, where Rb is 
the distance to the Galactic center ~ 8 kpc, and the 
observed fluxes from the disk are equal to the disk lu- 
minosities divided by only about 0.56x47ri?^ (Weiden- 
spointner et al. 2008a). Thus the broad line flux is 
(0.30 ± 0.08) X 10^3 ^ cm-2 g-^ from the bulge and 
(0.02 ±0.01) X 10"-'' 7 cm-2 g-i ^^^^ ^j^^ ^^^^^ ^j^^g ^-^^ 
narrow line flux is (0.44 ±0.08) x 10^^ 7 cm~^ s^^ from 
the bulge and (0.89 ± 0.14) x 10"^ 7 cm"^ s"^ from the 
disk. This gives a total observed broad 511 keV line flux 
of (0.32 ± 0.11) X 10""^ 7 cm~^ s~^ and a narrow line 
flux of (1.33 ± 0.16) X 10^'^ 7 cm^^ s^^, which gives an 
apparent broad/narrow line ratio of 0.24±0.08. This ra- 
tio would be observed, however, only if the full disk were 
observed, and for smaller solid angles the apparent ra- 
tio would be expected to increase approaching the bulge 
value of ^ 0.85, when only the bulge is viewed. Thus, 
we expect the apparent ratio to be quite sensitive to the 
effective sky coverage, and we suggest that this angular 
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dependence can test that prediction and further probe the 
nature and structure of the interstellar medium. 

Indeed, analyses of SPI/INTEGRAL observations of 
the Galactic bulge by Jean et al. (2006) give a broad 
line flux of (0.35±0.11) x 10~^ 7 cni~^ s~^ and a narrow 
line flux of (0.72 ± 0.12) x 10^^ 7 cm"^ s~\ which gives 
a broad/narrow line flux ratio of 0.49±0.17. A separate 
analysis by Churazov et al. (2005) gives a broad line 
flux of (0.24±0.10) X 10^^ 7 cm^^ s^^ and a narrow line 
flux of (0.51 ± 0.07) X 10^'^ 7 cni^^ s^^, which gives a 
broad/narrow line flux ratio of 0.47±0.20. Considering 
the large uncertainties in both the measurements and the 
model value, the expected bulge ratio of ^ 0.68±0.26 is 
well within 1 — tr. 

Moreover, within the positron bulge (< 1.5 kpc), we 
expect that the broad and narrow components of the 511 
keV line each essentially come from separate and poten- 
tially resolvable regions, and provide a further test. The 
narrow line emission is expected to come effectively all 
from central molecular zone and inner tilted disk within 
0.5 kpc, arising from thermal positronium and direct an- 
nihilation in the warm HII shells surrounding molecular 
clouds, ionized by FUV from central O stars. Essen- 
tially all of the broad line emission on the other hand 
is expected to come from the surrounding outer tilted 
disk between 0.5 and 1.5 kpc, arising from annihilation 
of positronium formed in flight in the cold neutral HI 
envelope of molecular clouds, where there is insufhcient 
FUV to generate surrounding HII shells. 

In particular, from the positron annihilation rates in 
Tables 2 & 6, we expect from the HII shells within 0.5 
kpc of the Galactic Center a narrow line flux of (0.39 ± 
0.06) X 10"'^ 7 cm~^ s^^ and no significant broad line 
flux, giving an expected broad/narrow line ratio of just 
^ within 0.5 kpc. Whereas from the HI gas in the 
surrounding tilted disk, we expect a broad line flux of 
(0.30 ± 0.10) X 10^'^ 7 cm~^ s~^, and a narrow line flux 
of (0.05 ± 0.03) X lO^'' 7 cm~^ s""'^ from direct thermal 
annihilation, giving an expected broad/narrow line ratio 
of ^ 6±3 between 0.5 kpc and 1.5 kpc, which is nearly 10 
times the combined full bulge average. Spectral analysis 
of the emission from these two regions should be able to 
test this prediction. 

7.6. Halo Contributions 

The 511 keV line flux from the halo makes up nearly 
half of the total Galactic annihilation line flux measured 
by SPI/INTEGRAL (Weidenspointner et al. 2008a) and 
from our estimates of positron propagation we also ex- 
pect that it accounts for a very similar (~ 49%) of the 
total line flux. Thus, as we discussed above, since this 
line intensity depends rather strongly on the mean dust 
content of the hot halo plasma, we expect that the Galac- 
tic averages of both the positronium fraction and the 
broad/narrow 511 keV line flux ratio are also very sen- 
sitive to the amount of dust in the halo. For as Jean et 
al. (2006) have shown and we list in Table 6, if the 
refractory elements in the hot plasma are all in dust 
grains, only ~ 18% of the positron annihilation occurs 
via positronium formation, and fully ^ 90% of the 511 
keV line emission is narrow with a FWHM ~ 2 keV, 
while the remaining 10% is very broad (FWHM ^11 
keV). This gives a halo flux of ~ (1.60 ± 0.37) x 10"^ 
cm~^ s~^ in the narrow line alone. That in turn gives a 



combined bulge and halo flux of ^ (2.34 ± 0.40) x 10~^ 
cm~^ s^^, which is quite consistent with that from the 
SPI/INTEGRAL best fit fiux of - (2.14 ± 0.11) x lO^^ 
cm~^ s~^ (Weidenspointner et al. 2008a), see Table 7. 
Combined with disk flux, that would give a total Galac- 
tic bulge, disk and halo flux of - (3.25 ± 0.42) x lO'^ 
cm~^ s~^ which is also quite consistent with the com- 
bined best-fit flux of ~ (2.87 ± 0.25) x lO^^ cm^^ g-^, 
from the SPI/INTEGRAL analyses. This also gives 
Galactic average values of the positronium fraction of 
only ^ 0.66 ± 0.12 and the broad/narrow 511 keV line 
flux ratio of only ^ 0.17 ± 0.04, for which the best-fit 
SPI/INTEGRAL values have not yet been published. 

If on the other hand, the grains have all disintegrated, 
the halo flux would still be almost as large, ~ (1.41 ± 
0.31) X 10^^ cm~^ s~^, but it would all be in the very wide 
(^ 11 keV) 511 keV line, so the Galactic average value 
of the broad/narrow 511 keV line flux ratio would jump 
to ~ 1.30 ± 0.29, and the positronium fraction would 
increase only slightly to ^ 0.74 ± 0.13. The actual dust 
content of the halo undoubtedly lies in between and since 
we expect the broad/narrow 511 keV line ratio from the 
halo beyond 1.5 kpc to be extremely sensitive, ranging 
from ^ 0.1 to 00, SPI/INTEGRAL measurements of that 
ratio should provide a unique measure of the halo dust. 

8. CONCLUSION 

In conclusion, we find that roughly half of the Galac- 
tic positrons from supernova generated radionuclei are 
produced in he bulge and inner disk within the stellar 
bulge of _R < 3 kpc, and half in the stellar disk beyond 
that. Essentially all of the bulge positrons are produced 
in the very hot plasma, through which they diffuse with 
a calculated mean free path of only 1% of the size of the 
region. Therefore, before they can either slow down and 
annihilate in the plasma or escape beyond, most of them 
are stopped and annihilate in the very dense (^100-1000 
H cm~^) warm outer shells of the molecular clouds that 
lie within i? < 1.5 kpc. The disk positrons, on the other 
hand, are produced about equally in the warm HI gas and 
the hot superbubbles, which blow out into the halo. We 
calculate that roughly half of the disk positrons, which 
either stream through the low density (~0.1 H cm~^) 
neutral HI at close to c, or diffuse through the hot super- 
bubble plasma, escape into the overlying halo before they 
can stop and annihilate in the low density (0.3-3 H cm~^) 
photoionized outer shells of clouds and superbubbles in 
the disk. Thus, positron propagation easily explains the 
larger annihilation flux from the bulge. Moreover, it ex- 
plains the high observed postronium fraction and narrow 
observed 511 keV line width, that both require annihila- 
tion predominantly in the warm ionized gas and not in 
the hot plasma in which nearly all of the positrons are 
produced, either in this or in other suggested sources. 

We expect these conclusions to be quite robust, since 
they are based on 1) extensive observations of galactic su- 
pernova rates and distributions, 2) observationally con- 
flrmed models of supernova yields of radionuclei and their 
decay-positron survival fractions, 3) wide ranging obser- 
vations of the properties and distributions of the various 
phases of the interstellar medium and magnetic fields, 
and 4) the widely used photodissociation region (PDR) 
models of molecular cloud shells. From these we have 
calculated the expected propagation of positrons in the 
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various phases, assuming streaming in the neutral phases 
and diffusion in the ionized gas and plasma phases, us- 
ing diffusion mean free paths in the latter phases calcu- 
lated from a two-component, phenomenological model of 
anisotropic turbulence, which is consistent with observa- 
tions of particle propagation in our best natural labora- 
tory, the interplanetary medium. 

Thus, as we show in Table 7, when positron propa- 
gation is considered, the positrons from the f3~^-deca.y 
chains of the radioactive nuclei, ^^Ni, ''^Ti, and ^^Al, 
produced in Galactic supernovae, can fully account for 



all of the features of the diffuse Galactic 511 keV 
and 3-7 continuum annihilation radiation observed by 
SPI/INTEGRAL. We have also predicted additional 
measurable features that can not only further test such 
an origin of the positrons but provide new information 
on the nature of the interstellar medium. 

This work was supported by NASA's International 
Gamma- Ray Astrophysics Laboratory Science Program, 
grant NNG05GE70G. 
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Fig. 1. — A schematic view of the Galactic Bulge, Disk and Halo assumed in these calculations, showing the stellar bulge and disk 
distributions of stars (+) and the assumed boundaries of the interstellar gas and plasma subdivisions, chosen for comparison with the 
SPI/INTEGRAL observational analyses. 
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^0.5 




Fig. 2. — The cummulativc fraction of total (thick sohd curve) Galactic positron production -F(< r) in the bulge and disk as a function 
of Galactic radius r from equations (3) and (4). The assumed distribution of positrons from decay of ^^Ni and ^''Ti, produced by SNIa and 
SNIp, are shown for the stellar disk (dot-dashed curve), bulge (thin solid curve) and a star burst near the Galactic Gentcr (dotted curve), 
while those from decay of ^''Al, produced by massive Wolf Rayet stars and core collapse supernovae are shown in the short dashed curve. 
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-l^K Plasma 







Fig. 3. — A schematic view of molecular clouds in the Central Molecular Zone of the Inner Bulge within 0.5 kpc of the Galactic Center, 
together with a cross section of the HII and HI envelopes surrounding the H^ core of a model cloud. 



Galactic Positron Annihilation 



33 



10 



-50° 




+50" 



■ Norma 

■ Scutum-Crux 

- Sagittarius-Carina 
Perseus 



-10 







10 



Fig. 4. — The Galactic spiral arms, determined from the study by Vallee (2005), showing the asymmetry, as viewed from our solar 
perspective, of the inner arms within 50° to either side of the positron bulge (±10°), measured by SPI/INTEGRAL.Thc asymmetry in the 
annihilation radiation from the Galactic disk is a product of the Norma arm dominating the negative longitudes due to its proximity to the 
inner region of the Galaxy, and the weakness of the Sagittarius-Carina arm due to its greater distance from the Galactic center from which 
the annihilation rate falls off exponentially from the Galactic center and its lower density/activity (Drimmel & Spergel 2001) compared to 
the other 3 arms. We find 511 keV line disk fluxes of (3.8±1.0)xl0~^ photons cm~^s~^ from the negative longitudes, and (2.3±0.6)xl0^^ 
photons cm~'^s~^ from the positive longitudes, and a ncgative-to-positive longitude flux ratio of 1.6it0.4. These values are quite consistent 
with the corresponding best fit (Weidcnspointner et al. 2008b) SPI/INTEGRAL line fluxes. 
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TABLE 1 
Assumed Properties of the Galactic Interstellar Bulge (<3.5 kpc) 







Inner Bu 


Ige 




IVI 


ddle Bulg 




Outer Bulge 


Galactic R 




<0.5 kp 









.5-1.5 kpc 




1.5-3.5 kpc 


Phase 


VH 


HII 


HI 


H2 


HM 


HI 


H2 


Bo 


Filling Factor / 


1 


0.0002 


0.0001 


0.0004 


1 


0.00004 


0.00005 


1 


Temperature T K 


(90 - 6) X lO'^ 


5000 


150 


30 


5 X 10^ 


20 


20 


1 X 10^ 


Density n H cm^^ 


0.04-0.007 


100 


1000 


1900 H2 


0.005 


1000 


1000 H2 


0.002 


Ionization Fraction a; 


1 


1 3 


X 10-^ 


10-5 


1 


3 X 10-*^ 


10-5 


1 


Slowing Down t^^ yr 


(1 - 6) X 10^ 


450 


90 


20 


9 X 10^ 


900 


150 


2 X 10^ 


Magnetic Field B fiG 


^100-17 


-.100 


-.100 


^100 


--12 


^12 


^12 


4 


Fluctuation 5B/B 


-.0.12-0.06 


1.3 


- 


- 


^0.03 


- 


- 


-.0.5 


Flux Tube l^ pc 


10-300 


^7 


-.5 


>7 


^1000 


^12 


~5 


1000 & 2000 


Max MHD Scale lo pc 


10-50 


- 


- 


- 


50 


- 


- 


75 


A II pc 


25-30 


0.0017 


- 


- 


45 


- 


- 


14 


Slowing Down I j pc 


2000-5000 


0.3 


12 


3 


7000 


12 


6 


7000 


iB/i^B +lsd) 


0.01-0.03 


0.96 


0.30 


>0.70 


0.07 


0.50 


0.95 


0.13 & 0.22 


Q^ 1043e+/. 


0.31±0.07 


^0 


-.0 


^0 


0.37±0.09 


^0 


^0 


0.56±0.14 
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TABLE 2 
Positron Propagation & Annihilation in the Galactic Interstellar Bulge (<3.5 kpc) 





Inner Bulge 


Ivliddle Bulge 


Outer Bulge 


Galactic R 


<0.5 kpc 


0.5-1.5 kpc 


1.5-3.5 kpc 


Phase X 


VH hi: 


HM HI 


Bo 


Filling Factor / 


1 -^0 


1 -^0 


1 


QX 10*3e+/B 


n.31±0.07 -^0 


0.37±0.09 -.0 


0.56±0.14 


Propagation 








Inner Bulge 








^VH:X 


^0.02 -.0.96 


- 


[-0.02] 


Middle Bulge 








^HM:X 


^0.22 


-.0.07 -.0.69 


[-0.02] 


Outer Bulge 








PBo:X 


^0.13 


-.0.04 -.0.40 


[- 0.14] 


Annihilation 








Inner Bulge 








Q^VH-.X 10'^3,+ /, 


^0.006 0.30±0.06 


- 


[ — 0.006] 


Middle Bulge 








Q^HM:X 


0.0S±0.04 


O.D3±0.01 0.26±0.13 


[-0.007] 


Outer Bulge 








Q^B.:X 


0.08±0.04 


0.02±0.01 0.22±0.11 


[ — 0.08] 


A^ 1043e+/. 


-.0.006 0.46±0.08 


0.05±0.02 0.4S±0.17 


[-0.09] 
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TABLE 3 
Assumed Properties of the Outer Bulge. Disk & Halo 





Outer Bulge 




Di 


sk 




Halo 


Galactic R 


1.5-3.5 kpc 




>3.5 


kpc 




>1.5 kpc 


Phase 


Bo 


HT 


HI 


WI 


HII 


H 


Filling Factor / 


1 


0.2 


0.5 


0.15 


0.15 


1 


Temperature T K 


1 X 10^ 


1 X 10^ 


8000 


8000 


8000 


1.6 X 10*^ 


Density tt H cm-^ 


0.002 


0.002 


0.13 


0.30 


2.6 


0.002 


Ionization Fraction x 


1 


1 


0.02 


0.9 


0.95 


1 


Slowing Down t^^ Myr 


23 


23 


0.7 


0.16 


0.018 


23 


Magnetic Field B fiG 


4 


4 


4 


4 


4 


4 


Fluctuation 5B/B 


-.0.5 


-.0.5 


- 


^0.5 


-.1 


-0.5 


Flux Tube l^ pc 


1000&2000 


<1000 


-.3000 


160 


85 


> 10000 


Max MHD Scale lo pc 


75 


75 


- 


10 


2 


75 


A|| pc 


14 


14 


- 


0.2 


0.1 


14 


Slowing Down l^^ pc 


7000 


7000 


90,000 
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15 


7000 


iB/i^B+^sd) 


0.13&0.22 


<0.13 


0.03 


0.56 


0.85 


>0.6 


QX 1043e+/s 


0.56±0.14 


0.45±0.12 


0.45±0.15 


0.13±0.04 


0.13±0.04 
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TABLE 4 
Positron Propagation & Annihilation in the Outer Bulge. Disk & Halo 





Outer Bulge 




D 


isk 




Halo 


Galactic R 


1.5-3.5 kpc 




>3. 


5 kpc 




> 1.5 kpc 


Phase X 


Bo 


HT 


HI 


WI 


HII 


H 


Filling Factor fj^ 


1 


0.2±0.05 


0.5±0.1 


0.15±0.05 


0.15±0.05 


1 


QX 10*3e+/B 


0.56±0.14 


0.45±0.15 


0.45±0.15 


0.13±0.04 


0.13±0.04 





Propagation 














^Bi/ B-m-.X 


[^0.04] 


- 


- 


- 


- 


— 0.04 


PBo:X 


[-0.14] 


- 


- 


- 


0.14±0.07 


0.29±0.10 


Pht.x 


- 


[<0.08] 


- 


0.15±0.07 


0.15±0.07 


0.70±0.35 


Phi.x 


- 


- 


— 0.02 


0.16±0.08 


0.16±0.08 


0.66±0.33 


PWI:X 


- 


-0 


— 


0.76±0.38 


0.05±0.02 


0.19±0.09 


Phii-.x 


— 


-0 


— 


0.03±0.01 


0.85±0.42 


0.10±0.05 


PH:X 


— n 


— 


— 


— 


— 


— 1 


Annihilation 














QPB^/Bm■.X lO^^e+Zs 


[-0.01] 


- 


- 


- 


- 


— 0.01 


QPbo-.x 


[—0.08] 


- 


- 


- 


0.0S±0.04 


0.16±0.08 


QPrt.x 


- 


[<0.04] 


- 


0.07±0.03 


0.07±0.03 


0.31±0.15 


QPhi.x 


- 


- 


-0.01 


0.07±0.03 


0.07±0.03 


0.30±0.15 


QPwi-.x 


- 


— 


— 


0.10±0.05 


O.OliO.Ol 


0.02±0.01 


QPrii.x 


— 


— 


— 


— 


0.11±0.05 


0.02±0.01 


QPH:X 


— 


-0 


— 


— 


-0 


-0 


A 10*3 e+ /a 


[-0.09] 


[-0.04] 


-0.01 


0.24±0.07 


0.34±0.0S 


0.82±0.23 
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TABLE 5 
Summary of Positron Production & Annihilation in Positron Bulge, Galactic Disk & Halo 





Production Rate Annihilation Rate 


Production 


Annihilation 




Q 10*3 e+/s 


A 10« e+/s 


Q/QTot % 


A/ATot % 


Bulge < 1.5 kpc 


0.68±0.11 


l.OOiO.19 


28±6 


41±8 


Bulge 1.5-3.5 kpc 


0.56±0.14 


[0.09]* 


23±6 


[4]* 


Disk > 3.5 kpc 


1.16±0.22 


0.59±0.11 


48±15 


25±5 


Halo > 1.5 kpc 





0.82±0.21 





34±9 


Hot Plasma H + 


1.66±0.24 


0.88±0.21 


69±13 


37±9 


Warm HII & Wl 


0.27±0.06 


1.04±0.14 


11±3 


43±6 


Warm Neutral HI 


0.45±0.15 


~0.01 


19±7 


~0 


Cold Neutral HI 


~0 


0.48±0.17 


^0 


20±7 


Molecular H2 


~0 


~0 


~0 


~0 


Hot Plasma 










Bulge < 1.5 kpc 


0.65±0.11 


0.06±0.02 


27±6 


3±1 


Bulge 1.5-3.5 kpc 


0.56±0.14 


[0.09]* 


23±6 


[4]* 


Disk > 3.5 kpc 


0.45±0.12 


[0.04]* 


19±5 


[2]* 


Halo > 1.5 kpc 





0.82±0.21 





34±9 


Warm HII & WI 










Bulge < 1.5 kpc 


~0.01 


0.46±0.08 


^0 


19±4 


Disk > 3.5 kpc 


0.26±0.06 


0.58±0.11 


11±3 


24±5 


Warm Neutral HI 










Disk > 3.5 kpc 


0.45±0.15 


~0.01 


19±7 


~0 


Cold Neutral HI 










Bulge < 1.5 kpc 


~0 


0.48±0.17 


^0 


20±7 
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TABLE 6 
Galactic Positronium (Ps) Fraction & 511 keV Line Luminosities & Fluxes 





Bi 


ilge 


Di 


3k 




Halo 




HII 


HI 


HII & WI 




HI 


Grains 


No Grains 




< 0.5 kpc 


0.5-1.5 kpc 


> 3.5 kpc 


> 


3.5 kpc 


> 1.5 kpc 


> 1.5 kpc 


Ps Fraction 


0.89 


0.95 


0.89 




0.99 


0.18 


0.42 


In-Flight Ps Fraction 


0.01 


0.95 


0.01 




0.70 








Thermal Ps Fraction 


0.88 


^0 


0.88 




0.29 


0.18 


0.42 


Direct Processes Fract 


ion 0.11 


0.05 


0.11 




0.01 


0.82 


0.58 


e+/T511 


1.50 


1.74 


1.50 




1.94 


0.58 


0.73 


Broad 7511 % 





83 







68 


10 


100 


Narrow 7511 % 


100 


17 


100 




32 


90 





Annihilation Rate 


10^3 e+ s-1 














Total Positron 


0.46±0.08 


0.48±0.17 


0.5S±0.11 




^0.01 


0.82±0.21 


0.82±0.21 


Total Positronium 


0.41±0.0S 


0.46±0.16 


0.52±0.10 




^0.01 


0.15±0.04 


0.34±0.09 


In-FHght Positronium 





0.46±0.16 







^0.01 








Thermal Positronium 


0.41±0.0S 


^0 


0.52±0.10 




-.0 


0.15±0.04 


0.34±0.09 


Direct Processes 


0.05±0.02 


0.02±0.01 


0.06±0.01 




-^0 


0.67±0.17 


0.48±0.12 



Mean Ps Fraction 
Luminosity 10"*' 
Ps 511 licV Line 
Direct 511 licV Lim 
Total 511 ItcV Line 
3-1 Continuum 
Broad 511 ItcV Lin< 
Narrow 511 ](eV Li] 
Broad 511 IteV Lin( 
Narrow 511 IseV Li] 
Broad/Na 
Broad/Na 



0.30±0.12 



0.20±0.04 
0.10±0.03 
0.30±0.05 
0.88±0.1S 

0.30±0.05 



0.23±0.08 
0.04±0.02 
0.27±0.08 
1.04±0.38 
0.23±0.08 
0.04±0.02 



^0 

~0.01 

0.02±0.01 

~0.01 

^0 



Ratic 
Ratic 



0.23±0.08 
0.34±0.05 
0.68±0.26 



Flu 



10" 



Total 511 ItcV Line 
3t Continuum 
Broad 511 licV Lin< 
Narrow 511 licV Li] 
Broad/Narrow Line 
Total 511 ItcV Line 
Broad 511 licV Lin< 
Narrow 511 IteV Li] 
Broad/Narrow Line 
Broad 511 IteV Lint 
Narrow 511 l^eV Li] 
Broad/Narrow Line 



0.39±0.06 
1.15±0.24 



0.39±0.06 
~0 



0.35±0.12 
1.36±0.50 
0.30±0.10 
0.05±0.03 
6±3 



0.74±0.14 
0.30±0.10 
0.44±0.08 
0.68±0.26 



0.26±0.05 
0.12±0.04 
0.38±0.06 
1.17±0.23 

0.38±0.06 

~0.01 

0.38±0.06 

~0.03 



0.89±0.14 ~0.02 
2.73±0.65 0.05±0.02 

~0 ~0.02 

0.89±0.14 ~0 

~0.02 
0.91±0.14 

~0.02 

0.89±0.13 

^0.02 



0.17±0.05 
0.96±0.24 
1.13±0.25 
0.77±0.19 
1.13±0.25 



0.07±0.n2 
1.35±0.34 
1.42±0.34 
0.33±0.08 
0.14±0.04 
1.28±0.30 

< 0.64 ± 0.15 > 

< 0.64 ± 0.15 > 
< 1.0 ± 0.2 > 



1.77±0.43 1.41±0.31 
0.41±0.11 0.45±0.24 
0.17±0.05 1.41±0.31 
1.59±0.37 

0.11 oo 

< 1.59 ± 0.37 > 

< 0.79 ± 0.19 > 

< 0.79 ± 0.19 > 
< 1.0 ± 0.2 > 



0.32±0.10 
1.33±0.15 
0.24±0.08 
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TABLE 7 
Expected & Observed Galactic Positron Annihilation Radiation* 



511 keV Line Flux 

10-3 ^./cm^ a 
Bulge Component 
Disk Component 
Halo & Bulge Componc 
Bulge Core <0.5 kpc 
Bulge Shell 0.5-1.5 kpc 
Bulge <1.5 kpc 
Disk -50° < I < 0° 
Disk +50° > ! > 0° 
Disk Ratio (-50° - 0° 



-0°) 



0.74±0.14 
0.91±0.14 
2.32±0.40 
0.39±0.10 
0.33±0.07 
0.72±0.13 
~0.3S 
^0.23 
-1.6 



0.75±0.09 
0.94±0.16 
2.14±0.11 
0.38±0.03 
0.41±0.06 
0.79±0.07 
0.43±0.05 
0.24±0.05 
1.8±0.4 



Wcidenspointncr ct al. 2008a 



Wcidcnspointncr et al. 2007 



Wcidenspointncr ct al. 2008b 



Annihilation Rate 

10^3 e+/s 

Bulge Component 

Disk Component 

Halo Sz Bulge Combined 

Bulge/Disk Component Rati< 

Positronium Fraction 

Inner Galaxy 



Disk > 3.5 kpc 
Halo > 1.5 kpc 

Broad & Narrow 511 keV Lii 

Bulge Warm Ionized Gas % 
Bulge Cold Neutral Gas % 
Broad/Narrow Line Ratio 
Bulge Region 

Bulge < 0.5 kpc 
Bulge 0.5-1.5 kpc 
Disk > 3.5 kpc 



Weidenspointner et al. 200Sa 



1.04±0.15 1.15±0.16 

0.76±0.11 0.81±0.14 

1.95±0.28 3.13[1.86*]±0.20 " [*adjusted c+ /-f^n] 



1.4±0.4 


1.4±0.3 


0.92±0.0S 


0.94±0.06 




0.95±0.03 




0.92±0.09 


0.88±0.08 


- ? - 


0.65±0.07 


- ? - 


49±12 


51±3 


51±14 


49±3 


0.68±0.26 


0.49±0.17 




0.47±0.20 


~0 


- ? - 


6±3 


- ? - 


~0.02 


- ? - 



Churazov et al. 2005 

Jean et al. 2006 

Weidenspointner et al. 2006 



I et al. 2006 



I et al. 2006 
'.ov et al. 2005 



For a Solar Galacto centric distance of 8.5 kpc, assumed in the INTEGRAL analyses, instead of 8 kpc, assumed in this study. 



